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We review the present state of understanding of solid friction at low velocities and for systems 
with negligibly small wear effects. 

We first analyze in detail the behavior of friction at interfaces between wacroscopic hard rough 
solids, whose main dynamical features are well described by the Rice-Ruina rate and state dependent 
constitutive law. We show that it results from two combined effects : (i) the threshold rheology 
of nanometer-thick junctions jammed under confinement into a soft glassy structure (ii) geometric 
aging, i.e. slow growth of the real arrea of contact via asperity creep interrupted by sliding. 

Closer analysis leads to identifying a second aging-rejuvenation process, at work within the junc- 
tions themselves. We compare the effects of structural aging at such multicontact, very highly 
confined, interfaces with those met under different confinement levels, namely boundary lubricated 
contacts and extended adhesive interfaces involving soft materials (hydrogels, elastomers). This 
leads us to propose a classification of frictional junctions in terms of the relative importance of 
jamming and adsoprtion-induced metastability. 



I. INTRODUCTION 



An established tradition, when writing about solid fric- 
tion, is to date its emergence as a well identified scientific 
question from Leonardo da Vinci (ca 1500). This makes 
it an unusually longstanding problem, since both its fun- 
damental physical aspects and its modelization for the 
purpose of studies of e.g. seismic fault dynamics are still 
under lively debate nowadays (see for example Q], 0, 
3). A very important progress in its phenomeno- 
logic modelling, for macroscopic solids, was accomplished 
with the formulation of the so-called rate and state con- 
stitutive laws, which emerged in the 70 's from the work 
on rock friction of Dieterich [j| , Rice and Ruina @ . 

It was later shown that, inspite of their simplicity in 
terms of number of parameters, they provide an excel- 
lent description of most of the salient features of the low- 
velocity frictional dynamics of a wide variety of materials, 
ranging from granite to paper. Such an amazing "uni- 
versality" is naturally appealing for the physicist, since 
it suggests : 

- The possibility of a unified, largely material- 
independent, description on the underlying microscopic 
level. 

- One step further, a possible feedback in terms of lim- 
itations and/or refinement of the mechanical constitutive 
laws. 

It is on this particular approach, deliberately different 
from that of tribology proper 1 , but which parallels those 
presently developed in the fields of plasticity of amor- 
phous materials |7j and rheology of jammed systems [8|, 
that we concentrate here. 



A. From Amontons-Coulomb to Rate-and-State 

For more than two centuries, the description of solid 
friction commonly used in mechanical modelling was pro- 
vided by the classical Amontons-Coulomb laws. They 
state that, when a nominally planar solid block lying on 
top of a planar track is submitted to a normal force W 
and a tangential one F (Figure ^1 : 

- No motion occurs as long as F is smaller than a finite 
threshold F s . 

- Sliding is dissipative, and the corresponding dynam- 
ical friction force Fd is constant and equal to F s . 

- Their common value F is proportional to the normal 
load W and, for a given W, independent of the macro- 
scopic contact area E. Hence, the frictional behavior of 
a couple of materials is characterized by a single number, 
the friction coefficient : 
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FIG. 1: A schematic slider-on-track system : the slider is 
pressed onto the track by the normal force W , and pulled 
by the tangential force F imposed through a driving stage of 
stiffness K . 



1 This means, in particular, that we limit ourselves to systems and 
sliding regimes such that wear and frictional self-heating play a 
negligible role. 



Note that this behavior (static threshold plus velocity- 
independent dynamic force) is the exact analogue, for 
friction, of the Hill rigid-plastic model of plasticity 
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(threshold plus constant yield stress), and thus proba- 
bly carries unphysical singularities as is the case of the 
Hill model. 

Indeed, various departures from this description or its 
implications have gradually emerged, the most salient of 
which are the following : 

(i) In general, the static friction coefficient fi a is larger 
than the dynamic one /id- 

(ii) n s is not a mere number, but a slowly increasing 
function of the so-called waiting time t w , i.e. the duration 
of static contact prior to sliding. 

(iii) When measured in stationary sliding, fid is not 
constant. In particular, for low enough velocities (typi- 
cally < 100/im.sec -1 ) it is a slowly decreasing function of 
V. 

That a stationary fid can be measured is then by itself 
a puzzle, since such a velocity-weakening characteristic 
is well known to make steady motion always unstable : a 
positive (negative) velocity fluctuation induces a decrease 
(increase) of the friction force, leading to further accela- 
ration (deceleration), thereby getting amplified. Indeed, 
when a slider is pulled at a low enough velocity through 
a driving stage, necessarily of finite stiffness, non steady 
sliding is frequently observed : the motion is jerky, al- 
ternating between "stick" periods of rest during which 
the stage stores elastic energy, and slip events (Figure 
EJ. However, as skilled mechanical engineers have known 
already for long, these stick-slip oscillations disappear 
upon stiffening the driving stage, and steady sliding is 
realized. 
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FIG. 3: Hysteretic friction force response of a paper/paper 
system in non-steady sliding (open circles) and for steady slid- 
ing at various velocities V (full circles). 



of which is controlled by a characteristic length Dq, of 
micrometric order. That is, its duration At « Do/Vf. 
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FIG. 2: Reduced tangential force vs time in the stick-slip 
regime for a paper/paper interface. From ref.|jj. 

This paradox points towards some inadequacy in the 
implicit assumption underlying the above instability ar- 
gument - namely that the stationary fid(V) function can 
be used as such to describe non steady motion. In other 
words, one suspects that the dynamic friction force does 
not only depends on the rate variable, i.e. the instanta- 
neous sliding velocity. 

(iv) Non steady friction is hysteretic, as illustrated on 
Figure [21 which shows the instantaneous friction force 
associated with a velocity cycle for which inertia is neg- 
ligible. 

(v) Dieterich p| has studied the frictional response to 
a sudden jump of the driving velocity from V$ to Vt (Fig- 
ure 0}. He showed that it exhibits a transient the span 



FIG. 4: Friction force transient following a jump of the driv- 
ing velocity at t = from Vi to Vf, for a PMMA/PMMA 
interface. (Courtesy of L. Bureau) 

These observations can be translated into the following 
statements : 

• A frictional interface at rest becomes stronger as time 
lapses : it strengthens when aging (see (ii)). 

• When sliding, it becomes weaker (see (i)) : it can be 
said to rejuvenate upon sliding. Moreover, its dynamic 
age in steady motion decreases with increasing V (see 
(iii)). 

• This, together with (iv) and (v), means that its phys- 
ical state evolves with a dynamics characterized by the 
length D and coupled to the sliding dynamics itself. 

Rate and state constitutive laws model such behaviors 
in terms of a single (or a few) dynamical state variable(s), 
the physical nature of which often remains unspecified. 
In such a situation, the physicist's approach aims at iden- 
tifying the physical content embedded in the state vari- 
ables, and at justifying on this basis their phenomenolog- 
ical dynamics. One may then hope, as a side benefit of 
importance, to be able to predict their limits of validity 
and, possibly, to propose some further extensions. 
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B. Spatial scales 



For this purpose, a natural first step is to identify 
the relevant length scales. In the case of friction 
between macroscopic solids, this immediately leads to 
distinguishing between two classes of systems. 

(A) Rough hard solids for which "reasonable" load- 
ing levels (apparent pressures W/H well below elastic 
moduli) do not result in an intimate molecular contact 
along the whole interfacial area. 

For them a first relevant length scale is the average 
size of the microcontacts - the spots which form the real 
area of contact S r . It lies usually in the micrometric 
range. Moreover, we will see that a detailed analysis 
leads to attributing ultimately frictional dissipation to 
elementary mechanical instabilities involving molecular 
rearrangements on the nanometric scale. 

(B) Soft and/or smooth solids which are able 
to get into intimate contact everywhere along their 
interface. For them, the nanometric scale retains its 
importance, while the mesoscopic one becomes irrelevant. 



In both cases, due to the long range of elastic inter- 
actions, a third important length scale is the global size 
of the system L. Its value has a crucial bearing on the 
spatial (in)homogeneity of the frictional motion, hence 
on the complexity of the sliding dynamics. 

An important consequence, not to be overlooked, re- 
sults straightforwardly from the identification of these 
spatial scales. Reducing the "state" of the interface to 
one or a few variables necessarily implies a statistical 
averaging, which can be meaningful only if performed 
on a large number of units. This entails that a rate 
and state phenomenology of friction is legitimate only on 
scales larger than a finite cutoff. In particular, when us- 
ing discretized (block) models for e.g. numerical studies 
of extended interface dynamics, the basic block size must 
remain much larger than either the average intercontact 
distance (class (A)) or at least the ultimate nanometric 
scale (class (B)). That is, we contend that the question 
of the regularity of the continuum limit of discretized 
models — a subject of debate in the field of mode II in- 
terfacial fracture along seismic faults [13 — though of 
course mathematically sound, is physically irrelevant : if 
it should turn out that a continuum description of the 
fracture head zone would imply lengths smaller than the 
intercontact distance, the interface could no longer be 
assumed to have homogeneous frictional properties, but 
should be explicily treated as a juxtaposition of frictional 
(contact) and non-frictional (non-contact) regions. 



C. Outline 

Following Tabor [J]), it is both useful and physically 
sound to express the friction force between two solids as 

F = a s .E r (2) 

where S r is the real contact area, and the stress a s is the 
so-called shear strength of the interface. 

Section II is devoted to friction at multicontact inter- 
faces between rough hard materials (MCI), an interfa- 
cial configuration which is prevalent when dealing with 
macroscopic bodies. We first analyze their geometry, 
then show that the geometric factor E r plays, for these 
systems, the part of a state-dependent variable, governed 
by what we will call the geometric age <f>. We then show 
that, for MCI, the rate dependence of the rheological 
factor a s can be assigned to local mechanical instabil- 
ities within the nanometer-thick molecular "junctions" 
forming the real contacts. On this basis, the Rice-Ruina 
constitutive law, which they originally formulated on a 
phenomenological basis, results as a good approximation 
for the low velocity frictional behavior of MCI. We sketch 
out its consequences in terms of the sliding dynamics of a 
driven spring-block system. We also point out its various 
limits of validity. 

These limits turn out to be of two different types. 

• On the one hand, within the physical framework 
which identifies "state" with geometric age, the phe- 
nomenological functional expression of £ r (</>) is neces- 
sarily aproximate, and its limits can be evaluated. 

• On the other hand, the detailed analysis of exper- 
iments brings to light the limits of the above physical 
framework itself. That is, it prompts the idea that a sec- 
ond underlying slow dynamics manifests itself through 
the rheological factor : geometric age is not the whole 
story, frictional contacts also have a structural age. 

This leads us to concentrate, in Section III, on the 
rheology of frictional contacts. We analyze it for various 
configurations (rough-on-flat MCI, surface force aparatus 
(SFA), extended soft contacts) corresponding to various 
confinement and compacity levels. In all cases, it ap- 
pears that structural aging/rejuvenation mechanisms are 
indeed at work in the contact-forming nanometer-thick 
interfacial junctions. The associated dynamics is quite 
strongly system-dependent. As a guide for future inves- 
tigation, we suggest a first level of classification which 
distinguishes between two main types of dissipative be- 
haviors, namely : 

(i) jammed junction plasticity, akin to that of soft 
glassy materials; 

(ii) adsorption-desorption controlled dynamics. 

MCI microcontacts belong to class (i), gel/glass con- 
tacts to class (ii). Both mechanisms probably con- 
tribute to friction in SFA contacts, while the case of elas- 
tomer/glass remains open to discussion. 
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II. MULTICONTACT INTERFACES 

The renewal of interest for solid friction in the past 20 
years, triggered by the work of rock mechanicians, was 
aimed at modelling friction along seismic faults. For this 
reason, it naturally focussed on MCI. It has resulted in 
bringing to light robust features shared by a wide variety 
of materials. We will see that this has led to the building 
a complete framework of description, namely : 

— a predictive constitutive law 

— an underlying physical interpretation opening onto 
further questions concerning contact rheology (Section 
III). 

This is why MCI friction is treated here extensively. 

A. Geometry of multicontact interfaces 

1. Surface roughness 



sized the importance of distinguishing between the ap- 
parent and the real areas of contact, T, app and S r , and 
of evaluating the average microcontact size a. It has 
emerged (see below) that typical values for E r /S app and 
a are respectively 10 -3 and a few microns. This means 
that a typical MCI is formed of a sparse set of micro- 
contacts with average separation ~ 100/im - a picture 
which has been confirmed by direct optical observation 
C fl^l and Figure 0. It is intuitively clear that increas- 
ing the normal load W results in decreasing the distance 
between the average surface planes (the so-called "clo- 
sure" ) , hence in increasing the number of microcontacts 
as well as the area of preexisting ones : S r increases with 
W, as illustrated by Dieterich's observations 0. Can 
one make the functional nature of this dependence ex- 
plicit? Does it depend on the statistical nature of the 
surface roughness and if so, how? Could it be that this 
dependence would explain the Amontons proportionality 
between F and Wl 



Consider two macroscopic solids (Figure QJ, referred 
to as slider and track, with nominally planar contact 
surfaces. Nominally means here that they are flat on 
large scales comparable with the macroscopic slider lat- 
eral size. It is well known that, due to the unavoidable 
presence of steps, atomic flatness cannot be realized over 
lengths much beyond the micrometric scale. An excep- 
tion is provided by lamellar solids, in particular cleaved 
mica for which this range may reach up to ~ 1 cm. Op- 
tically smooth surfaces, with r.m.s. roughness on the 
order of a few nanometers over 1 cm 2 , are uncommon, 
examples being provided by float glass and some highly 
polished metals. In general, natural or ground surfaces 
exhibit a r.m.s. roughness between 0.1 and a few microns 
on 1 cm 2 . 

Such random surfaces exhibit asperities with dis- 
tributed heights, so that physical contact between them 
occurs only at random spots, the microcontacts between 
load-bearing asperities. They form what we call multi- 
contact interfaces (MCI) (Figure [SJ. 
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FIG. 5: (a) Schematic view of a multicontact interface be- 
tween two nominally planar rough solids, (b) Optical vizual- 
ization of the light-transmitting real contacts at the MCI be- 
tween two transparent epoxy resin blocks. The contrast has 
been inverted. (Courtesy of O. Ronsin) 

It is Tabor [TlT |. in his pioneering work, who empha- 



2. The single microcontact 

As a first step, let us recall a few basic results from con- 
tact mechanics concerning a single contact. Let us focus 
on the a single pair of contacting asperities, modelled as 
elastic spherical caps with a common radius of curvature 
R and Young modulus E, pressed together by a normal 
force w acting along their intercenter axis. The problem 
of the resulting elastic contact was fully solved by Hertz 
and is equivalent to that between a rigid plane and an 
elastic sphere of radius R* = R/2 and Young modulus 
E* = E/2(l - i/ 2 ), with v the Poisson modulus (lH] . 
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FIG. 6: Hertz contact between an elastic sphere and a rigid 
plane. 

Its solution can be evaluated as follows. Due to the 
spherical geometry, the contact diameter a, the "com- 
pression" S (Figure EJ) and R* are related by a 2 ~ R*5. 
As can be checked from the full Hertz solution [l3| .the 
elastic energy is essentially stored within a depth on the 
order of the contact radius, hence the relevant strain level 
e ~ S/a, so that the average normal stress p ~ wja? ~ 
E*e r~j 5 /a. From this it immediately results that, di- 
mensionally 



5 



fwR*\ 1/3 r ( w 2 \ 1/3 fwE* 2 V 



The exact expressions (see 0]) only differ by multi- 
plicative constants of order unity. 

Now, most solids are truly elastic only up to a yield 
stress Y, beyond which they start deforming plastically. 
So, the Hertz expressions above lose validity when the 
maximum normal stress below the contact reaches this 
level. Upon increasing the normal load, the size of the 
plastified region increases, until it occupies a volume 
~ a 3 . At this stage the contact has become "fully plas- 
tic" and deforms so that the normal stress remains quasi- 
constant : p ~ H, where H is the hardness 2 of the softer 
material. At room temperature, for metals, the ratio 
H/E ranges in general between 10~ 2 and 10~ 3 , while 
for polymer glasses, H/E ~ 10 _1 - 10~ 2 . Imagine now, 
following Bowden and Tabor , that the apparent pres- 
sure Papp = WfE a pp is large enough for this regime to be 
reached in all the microcontacts forming a MCI. Then, 
the real area of contact is such that S r /S app ~ p app /H , 
and E r ~ W/H. 

So, in the fully plastic regime, E r is proportional to the 
normal load, and Amontons's law simply follows. Indeed 
the friction coefficient reads : 

F_ ^. v. 

W W H K ' 

Since, in this approximation, a s is the shear strength 
under the constant normal stress H, /i is effectively W- 
independent. However rough it is, this approximation is 
illuminating in several respects. 

On the one hand, the a priori surprising fact that the 
values of dry friction coefficients depend only weakly on 
the mechanical properties of materials and are commonly 
a fraction of unity can now be translated into the state- 
ment : interfacial shear strength values are roughly com- 
parable with bulk yield stress levels — a point which will 
be of qualitative importance later on. 

On the other hand, it permits to get a somewhat more 
precise idea about the minimum loading level necessary 
to form a MCI between given materials. Consider as an 
example a 1 cm thick steel plate under its own weight. 
Then p a pp ~ 10 3 Pa, while H ~ 10 9 Pa, hence S r /E app ~ 
10 -6 : a 1 cm 2 surface would form no more than about 10 
microcontacts of area 10/mi 2 (see below) - hardly enough 
to form a decent statistical set! 

Tabor's suggestion, which he developed for the case of 
metals, later gave rise to a number of discussions about 



2 On the basis of geometric arguments, one estimates that H ~ 
3Y 14]. 
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the relevance of the full plasticity assumption to wider 
classes of materials which also obey Amontons's law while 
being less ductile than metals. This opened onto the 
problem of modelling the elastic-plastic contact between 
random surfaces. 



3. Area of contact between random surfaces 

The crucial contribution to this question was made by 
Greenwood and Williamson [l5|. Their model reduces 
the characterization of each of the contacting random 
surfaces to : 

- the statistical distribution of asperity summit heights 
above some average plane cj>(z). These asperities are as- 
sumed to have spherical tips. 

- the asperity radius of curvature R, assumed to be 
unique. 

- the number of summits N on the apparent surface 

y 

^app ■ 

They show that the problem of contact between two 
such nominally flat elastic random surfaces can be cast 
into that of a "composite" random surface and a rigid 
plane [l3| . 




FIG. 7: Schematic representation of a Greenwood- 
Williamson interface (see text). 

Let d be the separation between the average plane and 
the rigid one (Figure |7J), n the number of microcontacts 
at this separation, assumed to be dilute enough to be 
treated as independent Hertz contacts. Each asperity 
with height z > d contributes to n : 

/>oo 

n = N dz<j)(z) (5) 

Jd 

The compression of a contact is (z — d), its area ■nR{z — d), 
and it bears the load (4/3)tt j R 1 / 2 (z - d) 3 / 2 (eq.@, 0) 
so that the real area of contact 



S r = iV/ dznR{z- d)<t>{z) (6) 

Jd 

and d is related to the given total normal load W by : 
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W = N dz (4/3)irR 1/2 (z - df' 2 4>{z) (7) 



The main physical content of the Greenwood - 
Williamson (GW) results already emerges under the 
schematic assumption of an exponential distribution 
4>o(z) — s _1 exp(— z/s) for z > and zero otherwise. 
One then gets : 



1 



W 



Hi ER 1 / 2 s 3 / 2 



R\ 1/2 W 



r: I — 



Both E r and the number of contacts are proportional 
to W , and individual contact sizes grow under increasing 
load in such a way that the average contact radius : 



1/2 



(9) 



remains a load-independent constant. 

GW have shown numerically that these properties are 
conserved, to a very good approximation, for a gaussian 
distribution of summit heights, in the load range, rang- 
ing over several decades, such that 1 < n < JV. A 
rough evaluation of N is obtained by representing Yi app 
as densely paved by average asperities of curvature 



and height s, so that N 
then translates into 



'Sapp/Rs. The upper limit 



E 



or S r <C E 



(10) 



Note that this condition also ensures that elastic interac- 
tions between microcontacts can reasonably be neglected. 

A considerable amount of work [l^ has been devoted 
to evaluating the two geometrical parameters involved in 
the GW model on the basis of surface topography mea- 
surements. Curvature determinations have been hotly 
debated, since they involve a second order derivative of 
the profile, which makes them strongly noise-sensitive. 
When measured with profilometers of lateral resolution of 
micrometric order, typical values of R for surfaces blasted 
or lapped with abrasive powders commonly lie in the 10 
- lOO^m range E3 EH E3> while r.m.s. roughnesses 
s ~ 1/im. Then the corresponding average contact radii 
a ~ 3 - 10/im, an estimate which has been confirmed 
from the Dieterich-Kilgore visualizations E2- 

Modern experimental improvements have revealed that 
surface profiles are in general more complex than was 
initially assumed by GW, exhibiting multiscale roughness 
distributions El- They are often modelled as sclf-affinc 



surfaces 3 . 

Persson |2jj has recently proposed a theoretical treat- 
ment of the resulting contact problem. His calculation of 
the contact area between a flat compliant medium and a 
rigid self-affine surface follows a renormalization scheme 
for the stress distribution for profiles at growing stages 
of magnification £. This function he proves to obey, for 
complete contact between the surfaces, a Fokker-Planck- 
like equation, with £ playing the part of time. The real 
case of incomplete contact is then simulated via a bound- 
ary condition eliminating from the area of contact the 
regions sustaining tensile stress. This enables him to ex- 
press the solution in terms of the elastic response of the 
compliant medium to imposed displacements specified all 
along its surface — while partial contact actually corre- 
sponds to a mixed situation of given displacements within 
contacts and given (null) normal stress elsewhere. It is 
therefore difficult to assess the influence of the ansatz on 
the main result of interest here, namely that the real 
contact area still obeys the Amontons proportionality 
S r ~ W. 

Note that Persson's result depends crucially on the 
physical existence of a lower cut-off of the spatial scales, 
which is in practice of mesoscopic order. Indeed, even 
in the highly improbable case where the fresh surfaces 
would remain fractal down to the atomic scale, the nano- 
asperities would get plastified, and thus smoothed out, 
under the very high pressures they would experience 
upon contact. 

This remark brings up the important question of the 
limit of validity of the pure elastic contact models. 
Clearly, when thinking in GW's terms, plastic deforma- 
tion of asperities starts coming into play when the aver- 
age pressure p becomes comparable with the yield stress 
Y. From equations © and © , this condition can be 
expressed in terms of the dimensionless plasticity index 



(11) 



which depends both on the topographic and the mechan- 
ical properties of the interface. 

The elastic regime is limited to the region ip < 1, for 
which GW have shown that the fraction of plastified con- 
tacts is negligible. Persson's more elaborate theory yields 
a comparable evaluation, with s and R in equation l|ll(l 
now understood to be those on the large scale (upper 
space cutoff). This last result takes into account consis- 
tently the above-mentioned plastification on small lateral 
scales. 

On the other hand, for ip larger than a few units, most 
microcontacts are in a state of full plastic deformation. 




3 The statistical properties of a self affine surface z(x, y) are in- 
variant under the scaling transformation x — > C,x,y — » C,y,z — + 
C^~ D i z, where Df is the fractal dimension. 
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Multiscale roughness loses relevance, and one may resort 
to the fully plastic version of the GW model. The matter 
which has flown being assumed to redistribute on a scale 
much larger than contact radii (no "piling up" ) , geometry 
imposes that the area of a contact involving an asperity 
with initial compression [z — d) is 2ttR(z — d). It bears 
the load ira 2 H, with H the hardness. Then, trivially 



V V i 

S r = _ a= V2Rs (12) 
H 

One may therefore assert the robustness of three char- 
acteristics of multicontact interfaces, namely, to a good 
approximation : 

- The real area of contact is proportional to the normal 
load, and independent of the apparent area T, app . 

- The average contact radius is load-independent, thus 
introducing a mesoscopic length scale a defined by the 
topography of the contacting surfaces. 

- In a wide variety of cases, a is of micrometric order. 

All these statements have been directly confirmed from 
the analysis, by Dieterich and Kilgore 12], of their optical 
images of several MCI between transparent solids. 

For metals, E/Y ~ 10 2 — 10 3 . It would therefore take 
unrealistically small values of s/R, of order 10 -4 — 10~ 6 , 
for metal/metal MCI to be in the pure elastic regime. 
They are fully plastic, as anticipated by Bowden and Ta- 
bor. 

The opposite case is that of elastomers, which are fully 
(visco)elastic up to strains > 1. Their interfaces, such as 
that between tire and r oad, when analyzed in the frame 
of Persson's predictions [2(| |2l| should provide a test of 
his theory. 

MCI involving polymeric glasses {E/Y ~ 10-100) 
pertain to the intermediate regime, where ip > 1 : a 
fraction only of the microcontacts flow significantly. 
Although no quantitative theory is available, we believe 
that, in view of their robustness, the above-mentioned 
three main results of the GW model hold for all non- 
purely elastic MCI. Confirmation of this statement has 
been provided by experimental investigations of various 
interfacial properties, such as shear stiffness, which all 
exhibit the Amontons linear dependence on normal load 
(see Appendix A). 

In fine, we can summarize the above results into the 
following simple statement for the friction coefficient of 
a MCI : 



p is the average normal stress borne by the microcon- 
tacts. It is load-independent. Its expression in terms of 
the material properties (Young modulus and yield stress) 
and of the topographic surface characteristics depends 
upon the value of the plasticity index ip (eq.JTIJ). For 



ij) 3> 1 (fully plastic regime) p — H ~ 2>Y. In the oppo- 
site, fully elastic case ip -C 1, p ~ EyJ s/R ~ H/tp, and 
in the intermediate elasto-plastic range it is expected to 
extrapolate smoothly with ip between these two limits. 

Let us insist that the interfacial shear strength a s in 
expression ifHfl) . which is likely to depend upon the con- 
tact pressure, is load-independent for a given MCI, since 
such is p itself. This result is specific of multicontact in- 
terfaces, for which Amontons's law is of purely geometric 
origin. It must be contrasted with the case of intimate 
single contacts such as those studied with the surface 
force apparatus (SFA) or extended gel glass ones (see 
Section III). For these systems, Amontons's law, when 
observed, must be attributed to the pressure dependence 
of cr s itself [13. 

B. Geometric age : a major state variable for MCI 

1. Time dependence of the static threshold: 

The classical description of solid friction states that 
there exists, for any given interface, a well defined static 
friction coefficient fi s — F s /W such that, as long as the 
applied shear force F < F s , no sliding occurs. 

This amounts to asserting that, at a depinning thresh- 
old F s , the interface commutes from a purely elastic, re- 
versible response to external shearing, to an irreversible, 
dissipatively flowing one. Most likely, and in view of un- 
avoidable disorder and temperature effects, such a strong 
statement, which implies that the interface would un- 
dergo an abrupt unjamming transition, is only approxi- 
mate. We come back to this point in § II. C. 4, where we 
discuss in detail the intrinsic difficulties and limits asso- 
ciated with the definition of such a threshold. Let us only 
state at this point that, in order to minimize ambiguities 
and make comparison between data meaningful, it is im- 
portant to define carefully the protocol used to measure 
Ms- 

A series of such stop and go experiments consists in : 

(i) Preparing the initial interfacial state reproducibly, 
by sliding steadily at a chosen velocity V prep , then stop- 
ping the pulling. 

(ii) Waiting at rest for a given waiting time t w under 
a specified shear stress. This may either be self-selected 
by the system (natural arrest stress) or imposed at some 
value below this level. 

(iii) Resuming loading at a prescribed velocity Vi oa d, 
for example Vi oad = V prep . 

A typical shear force response is displayed on Figure|Sl 
fi s is then conventionally identified with the peak level. 
Experiments of this type have been perfored on a num- 
ber of MCI. They reveal that fi s is not a mere number 
characterizing a couple of solids, but a slowly increasing 
function of the waiting time. More specifically, fi s (t w ) 
varies logarithmically over several decades of t w rang- 
ing up from about 1 sec (the typical fast limit for such 
mechanical experiments). This behavior, displayed on 
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FIG. 8: Friction trace (reduced friction force /j, = F/W vs 
time t) in a stop and go experiment (see text). V pr e P ~ Vi oa d- 
Interface : PMM A/PMM A at T = 300K. 



Figure holds for a wide variety of materials including 
metals |23l , rocks p| jl^l , glassy polymers and paper 
0. 

As a MCI ages at rest, it strengthens logarithmically. 
Moreover, the slope 



D 



d(ln t v 



(14) 



is found, for all the above-mentioned materials, to be 
roughly (see Figure^!) on the order of 1CP 2 . 

Such a "generic" behavior is striking, and suggests that 
it results from a robust physical mechanism. 



comparable with the bulk yield stresses of the contacting 
materials 4 . At this stress level, which prevails in a vol- 
ume ~ a 3 (with a the microcontact radius), one expects 
materials to undergo plastic creep, resulting in the slow 
growth of microcontacts, hence of S r . 

Confirmation of this idea was obtained by Dieterich 
and Kilgore 01 wn0 were able to measure directly the 
time evolution of XV on their optical images of the mi- 
crocontacts forming the MCI between two transparent 
glassy polymer blocks. They found that £ r (i w ) docs 
grow logarithmically, at a rate compatible with that mea- 
sured in microindcntation experiments. 

The question is then to ascertain whether or not this 
area growth is sufficient to explain that of fi s . This 
has been investigated in detail by Berthoud et al |25j . 
They studied the temperature dependence of the static 
aging slope B for symmetric MCI involving the polymeric 
glasses PMMA and PS (polystyrene) between room tem- 
perature and the vicinity of the bulk glass transitions. 
They analyzed their results in terms of a model for the 
growth of E r due to Brechet and Estrin [2^ . 

This model schematizes the creep-induced growth of a 
microcontact as follows. Once a microcontact has been 
created, at time t = 0, in a first stage of fast plastic flow, 
of duration negligible on the scale of the later evolution, 
the normal stress a sets to an "initial" value oy of the 
order of the yield stress in the relevant geometry 5 . Once 
this state is reached, plastic evolution continues via creep, 
the rate of which is given by an expression a la Nabarro- 
Hcrring : 



SANDSTONE 





FIG. 9: Static friction coefficient fi s vs waiting time t w for 
symmetric MCIs. Data from : Heslot et al @ (paper), Baum- 
berger [2(J (PMMA), Dieterich ^ (sandstone), Dokos ^ 
(steel) . 



e = e exp a/S 



(15) 



where e is the compressive strain rate (treated as a 
scalar), S the so-called strain-rate sensitivity of the flow 
stress, and e a cr-independent Arrhenius factor. Since 
this creep law reflects a thermally activated process, both 
eo and S are T-dependent. 

Such plastic deformation occurs at constant volume, 
so that : 



1 d(a 2 ) 
a?, alt 



(16) 



with a the contact radius, do its initial value. Since 
a = w/ira 2 ~ (1 — €)w/ira^ + C(e 2 ), with w the normal 
load. Then, from eas. (|14fl and (|15fl . and approximating 
the MCI by a set of identical average GW contacts, the 
real area of contact evolves as : 



2. Creep growth of the real area of contact 

Tabor's decomposition of the friction force (eq.JSJ), 
together with the Greenwood- Williamson analysis, im- 
mediately raises the possibility that this strengthening 
might be attributable to slow growth of the real area 
of contact. Indeed, we have seen that, in general, the 
average normal stress p borne by the microcontacts is 



I + to In ( 1 + - 



G(ln 2 ) 



(17) 



4 More precisely, for an asymmetric interface, p is on the order of 
the yield stress of the softer material. 

5 In uniaxial loading, ay = Y, while for a sphere-sphere contact 
ay = H ~ 3Y 
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with 



m = — , r = e — e 1/m (18) 
ay rn 

Typically, at room temperature, for the systems of in- 
terest here, m values are on the order of a few 1CP 2 . On 
the other hand, no reliable evaluation of the cross-over 
time r can be made a priori, due not only to its exponen- 
tial dependence on m but even more to the lack of any 
precise data on eq. All MCI studied up to now, except for 
polymer glass ones close to the glass temperature, exhibit 
a linear logarithmic growth of /i s , hence of E r , for wait- 
ing times above 1 sec. This only enables us to state, at 
this stage, that in general r is smaller than this limiting 
value. Then, assuming that the Tabor interfacial shear 
strength a s exhibits no aging, one obtains from ea. (|17|l . 
in the accessible range t ^> r, for the log-slope of the 
static friction coefficient : 

B = /i s0 m (19) 

where /j s o = a s T, r o/W = o~ s /o~y is the static threshold 
base value at short times t <§; r 6 . Since friction coef- 
ficients of MCI are usually a fraction of unity, one thus 
expects, at room temperature, B to lie in the 10~ 2 range. 

The analysis of experimental data for B(T) in ref.[25j 
is quite intricate, in view of several difficulties concerned 
with : 

- the fact that r and, hence, fj, s o , could be accessed in 
the realizable waiting time range only close to the glass 
transition temperatures T g . 

- the problems related with mapping the authors' bulk 
data for uq and S onto the sphere-sphere geometry. 

Berthoud et al were nevertheless able to show that 
expression (tT§|) accounts semiquantitatively for the 
roughly tenfold increase of B between its values, of order 
10~ 2 at room temperature as expected from expression 
(|19fl . and the vicinity of the bulk T ff 's. However, it 
must be noted that ea. lfn?|l is found to systematically 
underestimate the experimental data, the relative misfit 
increasing significantly on approaching T g . This hints 
towards the fact that, while creep- induced growth of the 
real area of contact is responsible for most of the static 
strengthening, some aging of the interfacial strength is 
not completely ruled out : for example, close to T g , such 
symmetric polymer glass contacts might exhibit partial 
"healing" due to interdiffusion of polymer chains. 

Clearly, the Brechet-Estrin model only describes plas- 
tic creep of the contacting asperities. However, quasi- 
logarithmic static aging was also observed by Ronsin ct al 
[30| on rough rubber/rough glass MCIs. For such materi- 
als, contact area growth is obviously of viscoelastic origin. 



Strictly speaking the value of a s which appears in this expression 
should be that at the driving velocity (see §11. C. 3 and 121-1 



We show in Appendix B, following Hui et al j^l), that the 
Greenwood- Williamson model can be worked out explic- 
itly for linear viscoelastic materials, with a very simple 
outcome. Namely, the GW results (eq.©) for the elastic 
MCI still hold formally, provided that the inverse Young 
modulus l/E is replaced by the so-called creep compli- 
ance J(t), which measures the delayed strain response to 
a unit instantaneous stress jump. Hence the real area of 
contact becomes a slowly increasing function of contact 
duration. A logarithmic increase is known to provide a 
reasonable approximation for J{t) over time decades, for 
materials with a very wide spectrum of relaxation times 
such as rubbers. 



3. Geometric age as a dynamical state variable 

In the light of the foregoing analysis we will for the mo- 
ment attribute the strengthening of MCI static thresh- 
olds to the sole time dependence of the area factor in 
Tabor's expression. S r thus becomes a function of the 
geometric age of the interface (f> by which t should now 
be replaced in equation (|17|l . More precisely, this age is 
defined as follows. 

• For a non-moving MCI, <p(t) is simply the time which 
has been spent at rest at time t. 

• Consider now a MCI sliding at the constant veloc- 
ity V. As motion proceeds, a given microcontact, once 
created (born) is gradually sheared until it slides, then 
disappears (dies) when the relative displacement between 
the partner asperities reaches a value on the order of a 
fraction of the contact diameter. Since, under constant 
normal load, the average number of microcontacts is con- 
served, any contact death is, on average, associated with 
the birth of of a new microcontact at an uncorrelated po- 
sition. Since E r <C S app , we can safely consider that the 
newborn is formed between fresh asperities, which have 
not yet experienced creep 7 . 

Contact renewal therefore limits the age of the MCI to 
the average lifetime of a given configuration of microcon- 
tacts, which can be written phenomenologically as : 

4>ss = ^ (20) 

In other words, motion interrupts aging, since interfa- 
cial configuration memory is destroyed after sliding the 
characteristic length Dq, which we expect to lie in the 
micrometer range. 

Note that the larger V, the younger the steady sliding 
MCI, hence the smaller £ r (^>): geometric aging thus im- 
mediately appears as a candidate process for explaining 



This approximation, which neglects wear, is validated by the 
observed stability of the frictional characteristics of a MCI over 
slid lengths as large as a few 10 cm |32| — a distance over which 
we can estimate the number of contact configuration renewals to 
lO 4 " 5 . 
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the V-weakening behavior of the steady sliding dynami- 
cal friction coefficient Hd(V) mentioned in Section I. 

• Let us now turn to non-steady sliding at the instan- 
taneous velocity x(t). 4> is no longer time-independent, 
since it keeps track of geometric aging over the time nec- 
essary to slide the memory length Dq. In this sense, -Do 
is the length over which memory of the history of motion 
is preserved. 

The most simple phenomenological expression for <f> 
accounting for this behavior is : 



4>(t) = / dti exp 



(x(t) - x(ti) 



D 



(21) 



where to is the time at which the two solids were first 
brought into contact. Note that, in agreement, as needed, 
with the previously defined expressions. : 

- for the static MCI equation (|21[1 yields cj)(t) — (t — 
to) = tw, 

- steady sliding corresponds to the limit to — > — 00 , so 
that 4>(t) as defined from ea. <|21(l reduces to the constant 
<f> 8S (V) =D Q /V, 

Geometric age thus becomes a dynamical variable, cou- 
pled to the instantaneous velocity x(t) by the non-linear 
differential equation, equivalent to expression l(2*TJl: 



= 1 X(j) 
Do 



(22) 



which was first proposed by Rice and Ruina . 



That the memory of the interfacial state of sliding 
MCIs is indeed characterized by a length, a possibility 
first suggested by Rabinowicz [33j in an often overlooked 
pioneer work, was established by Dieterich [£j on the ba- 
sis of his systematic exploration of frictional transients 
following velocity jumps. In these experiments the slider 
is set into steady motion by driving it, through a "spring" 
of stiffness K 8 , at an initial velocity V = Vl. At t = 0, 
V is suddenly jumped to Vf, and one measures the force 
response. As can be seen on Figure 0] this exhibits a 
two-step transient. 

(i) In a first, very rapid stage, for upward (resp. down- 
ward) velocity jumps, the instantaneous dynamic friction 
coefficient increases (resp. decreases). 

(ii) This so-called "direct effect" is followed by a much 
slower monotonous variation in the reverse direction, 
ending at the level fid(Vf). 

Dieterich found that the duration of the transient, 
which is dominated by that of stage (ii), scales as 
l/Vf. From this he identified a value of Do — 5/xm 
for a granite/granite MCI. One may then conclude 
that the slow part of the transient corresponds to the 
gradual relaxation of the real contact area from its 



initial steady value E r (</> = Do/Vi) towards its final one 
S r (^ = D /F / ). 

In conclusion of this analysis, it appears reasonable to 
write tentatively Tabor's expression as : 



F = S r {<j>) a s (x) 



(23) 



which assigns the whole state dependence to the area 
factor, while assuming that the shear strength a s only 
depends on the instantaneous sliding velocity. The valid- 
ity of this assumption will be discussed at length in the 
forthcoming sections. 

Equations (|23|l and (|22H . together with Newton's equa- 
tion : 



M'x = Fdrive {X,t) — F (</>, i) 



(24) 



where Fdrive is the externally applied driving force, pro- 
vide a closed set of equations for the frictional motion of 
the slider once the form of the functional i-dependence 
of F has been identified. 



C. Junction rheology : gross features 

1. Junctions at multicontact interfaces 

Up to now, we have concentrated on the analysis of the 
contact geometry of MCI which, though an important 
prerequisite, does not yet touch upon our main question, 
namely that of the origin of frictional dissipation and of 
the detailed nature of the associated rheology described 
by the interfacial shear strength a s . Let us first try to 
identify the regions which are the seat of frictional dissi- 
pation. 

We have restricted our definition of MCI to the normal 
load range such that microcontacts are numerous enough 
to form a good statistical set, but sparse enough for elas- 
tic interactions between them via the bulk materials to 
be negligible. This entails that we can now simply focus 
on the behavior of a single sheared microcontact. Such 
a unit (Figure 110(1 is constituted of the bulk of the two 
asperities and of an interfacial layer in which molecules 
from both surfaces have come into adhesive contact. 




FIG. 10: Each microcontact is formed by a nm-thick adhesive 
junction between the two asperity bulks. 



More precisely, K is the equivalent stiffness of the driving stage 
plus slider system. 



This layer, which we call the junction, has a disor- 
dered structure. This is obvious when the two solids are 
polymer glasses, since in this case the junction is formed 
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by polymer tails and loops protruding from the amor- 
phous bulks. In the case of crystalline bulks, structural 
and chemical surface disorder, which prevails under usual 
(non atomically planar and non ultra high vacuum) con- 
ditions together with plastic smoothing out of nano-scale 
roughness certainly lead to structural disorder. 

Such highly defective structures, probably less dense 
than the bulks, are likely to result in lower resistance 
against plastic shear deformation. So, we will assume 
that the junction is a disordered quasi bidimensional 
medium with thickness h in the nanometric range, where 
shear naturally localizes. We will show that this assump- 
tion is borne out by the analysis of experimental data on 
low velocity friction. 



2. A threshold rheology 

Stating the existence of a well-defined static friction co- 
efficient can be formulated equivalcntly as the fact that, 
when the sliding velocity V — > 0, a s does not vanish, but 
tends towards a finite limit a*. That is, from a rheolog- 
ical viewpoint, frictional contacts behave as yield stress 
fluids. 

The yield stress a* therefore appears as the threshold 
beyond which the interfacial quasi-2D solid flows plasti- 
cally. 

At lower stress levels, in this simple picture, the inter- 
face does not slide, it is pinned and responds to shear 
as a solid i.e., in principle, clastically. Due to the nano- 
metric thickness of the junctions, the corresponding stiff- 
ness is always much larger than that of the microcontact- 
forming asperities themselves (see Appendix A) 9 , it is 
therefore not accessible in MCI configurations. 

It is common knowledge that plastic flow does not set 
in as an ideally sharp transition : at non-zero tempera- 
tures, thermal activation induces creep below the nomi- 
nal yield stress, which can be attained, in principle, only 
by loading at extremely high rates. In other words, the 
smaller the loading rate, the more fuzzy the threshold. 
So, one expects that a constant shear load close below 
the static threshold should induce creep-like frictional 
sliding. This was indeed observed by Heslot et al 0- 
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3. Beyond threshold : rate effects 

a. Velocity jumps : the direct effect. We already 
mentioned in §11. B. 3 that the force response (see Fig- 
ure 0} to a jump of the driving velocity from Vi to Vf 



9 A rough estimate for the corresponding ratio is a/h > 10 3 , with 

a the average microcontact radius, h the junction thickness. 
10 Note, however, that with MCI this creep is certainly amplified by 
the geometric rejuvenation (weakening) associated with sliding. 



exhibits, previous to the slow transient attributable to 
geometric age adaptation, a first much faster part. The 
associated force jump WA/id is positive (resp. negative) 
for Vf > Vi (resp. Vf < Vi), hence the term "direct 
effect" . Dieterich has shown that : 



AfidiVi ^ V f ) = A]n -f 



(25) 



where A is a constant for a given couple of materials. Its 
values lie in the 10 -2 range. This holds for all the MCI 
studied up to now, for velocities between about 0.1 and 
lOO/im.sec -1 . 

Note that, as long as inertia is negligible, motion is 
quasi-static, i.e. after the jump at t = : 



K(V f t-x) = W[ f , d (<t>,x)- f i s d t (V i )] (26) 

x(t) measures the position of some reference point 
on the slider along the pulling direction, x(0) = 0, 
and H d l iy) = [id (Dq/V, V) is the friction coefficient for 
steady sliding at velocity V. 

The peak value of the transient force (dfid/dt = 0) 
therefore occurs for x = Vf. Velocity jump experiments 
are performed with as stiff as possible driving stages such 
that, as can be checked by direct measurements, the dis- 
tance slept during the fast part of the transient be much 
smaller than the memory length Dq. Under such condi- 
tions the corresponding variation of geometric age from 
its initial value 0, = Dq/Vi is negligible, and the direct 
effect is fully attributable to the rate dependence of the 
shear strength a s . 

From equation (|23|l one then expects that : 



A/ x (Vi - V f ) = Er( y^ } [a 8 (V f ) - a.m (27) 
That is, making use of equations (jl7(l and l|25(l 



Aa s = a s (V f ) - a s (Vi) = 



A\n(V f /V) 



I ■!■ /// In ( 1 + ^ 



Since m ~ 10 2 , the log term in the denominator can 
be neglected, and Aa s ~ In (¥f:J = -Aay io- (^-J , 
from which we can write the following empirical expres- 
sion for <7, : 



a s (x) 



0~sQ 



1 + a In ■ 



V 



G(ln 2 ) 



(29) 



with 



Ao Y 

&s0 



(30) 
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and a s o = <t s (Vo) is the shear strength at the reference 
velocity Vg, which may be chosen anywhere in the range 
(0.1 - 100/zm/sec) where equation (|2*5|l holds. 

The necessity of resorting to an expression involving 
a finite reference velocity is imposed by the formal di- 
vergence of the logarithm in the vanishing x limit. This 
divergence is of course unphysical and only points out the 
limits of the empirical approach. We will return to this 
point later, when discussing the behavior to be expected 
for o~ s (x) in the very low and large velocity regimes. 

Note finally that (t s q/(Jy is the friction coefficient fi s o 
appearing in equation (|19|) , it is therefore roughly of order 
1 and, in order of magnitude : 

a-A-KT 2 (31) 

b. Steady sliding friction coefficient. The velocity 
jump experiments initiated by Dieterich have played a 
pioneering role to evidence the state and rate character 
of MCI friction and to separate clearly the age (state) and 
rheologic (rate) contributions to the dynamic friction co- 
efficient fid- However, it is the studies of the velocity and 
temperature dependences of fid(V) under steady sliding 
conditions which have led to a systematic confirmation 
of the validity of expression l|29[l and opened the way to 
its more quantitative explicitation. 

The typical experimental Independence of fid is illus- 
trated on Figure 1111 on the cases of granite, paper and 
PMMA symetric MCI at room temperature. It is seen 
that : 

(i) fid is V- weakening (dfid/dV < 0) in the whole low 
velocity regime explored in such experiments. 

(ii) Its variations are quasi-logarithmic, except in the 
higher V range (V ~ 100/^m/sec) where it exhibits a 
saturating trend. For all materials, the log-slope (3d — 
—dfidl 'd(lnV) lies, once more, in the 10 -2 range. 

Such results can be confronted with the predictions 
obtained from equations Ij23|l . (|17|l and (|29(l which yield, 
in steady motion where the geometric age <fi = Dq/V : 



i + m m(i + ^) 


. v 


+ a m — 




V 



where we have taken advantage of the smallness of m and 
a to neglect the C(ln 2 ) terms, and made use of equations 
(O and PJ|). 

For V <C Dq/t, equation lj3~2"|) gives for the log-slope of 
fid the value : 

p d = B-A (33) 

The fact that (3d is positive entails that, at room tem- 
perature, B > A. Moreover, the (3d data are found to 
be quantitatively compatible with the available indepen- 
dently measured B and A. 
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FIG. 11: Dynamic friction coefficient fid versus steady slid- 
ing velocity V for symmetric MCI : (a) granite (data from 
0) (b) paper(data from @) (c) PMMA at T = 296° K (full 
symbols) and T = 384° K ~ T g (empty symbols) (data from 
[23). Due to the increase of the crossover time r (see eg, 1181 
with T, close below T g , V m in has decreased enough to lie in 
the experimental range. 



On the other hand, the saturating trend at the larger 
V's appears as resulting from that of the real area of 
contact in the small age limit, which becomes relevant 
for V > Dq/t. More precisely, equation (1.'S2[| predicts a 
broad minimum of fid, which should occur for V m i n — 
(D /t)(B-A)/A^D /t. 

That is, expression ll-il't predicts that, beyond V m i n , 
dynamic friction should become velocity-strengthening. 
This natural physical consequence of the nature of ge- 
ometric aging, which has been checked experimentally 
on the case of paper , has been overlooked up to now 
in most studies of fault dynamics. It is important in view 
of its bearing upon the nature of the sliding dynamics. 

We have seen in §11. B. 2 that parameter B varies with 
temperature in a way which is explained by its physi- 
cal origin - the thermally activated creep of S r . On the 
same polymer glass systems, Berthoud et al have shown 
|29| that A increases quasi-linearly 11 when increasing T 
up to the vicinity of T g . Blanpied et al [35j and Nakatani 
|36j | have investigated the behavior of, respectively, (3 d 



1 Note, however, that the T-range in these experiments is small. 
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and A for granite over wide T-ranges, from room tem- 
perature to 800°C (still well below the melting point). 
Pd was found to decrease with increasing T, enough for 
hot granite to become velocity-strengthening - a point of 
importance in the context of deep seismicity. A, on the 
contrary, increases in a quasi-linear manner. We defer a 
more detailed analysis and interpretation of these results, 
in particular those concerned with A, to § II. C. 4 below, 
where we propose a physical model for the rate effect. 



4- Threshold rheology as a signature of multistability 

The primary feature of frictional rheology is the exis- 
tence of a threshold. That is, the force F needed for an 
interface to slide at velocity V remains finite for van- 
ishing Vs. This is obviously at odds with the stan- 
dard picture of dissipation at low rates, namely : if the 
externally imposed rate of shear is much smaller than 
the internal relaxation ones, the system evolves quasi- 
adiabatically, so that dissipation vanishes linearly with v. 
As first formulated by M. Brillouin in 1904 ^5%, in order 
for F to exhibit a finite threshold, it is therefore neces- 
sary that, however slow the drive, the sheared medium 
evolves through a succession of adiabatic adaptation pe- 
riods interspeded with fast instability events. Each of 
these events then corresponds, even at vanishing v, to a 
finite energy loss, thus accounting for the threshold be- 
havior. 

In modern language, a threshold rheology implies 
multistability, as has been developed at length in various 
fields such as magnetic hysteresis, wetting dynamics, 
char ge d ensity wave and type-II superconducting trans- 
port |38j . 

a. A toy model for junction rheology 

Let us first illustrate this idea here on a toy model of 
a sheared junction 12 . We represent it, as shown on Fig- 
ure El as formed by a set of identical blobs of elastic 
matter attached to two stiff plates with vertical separa- 
tion h, and area E. The blobs are randomly distributed 
along the surfaces of these plates, and form n contacts 
per unit area. They compress each other when making 
contact, in which case they interact via the (repulsive) 
pinning potential U(x), where x measures their relative 
position along the shearing direction (Figure O. As x is 
increased, some contacts are being destroyed while oth- 
ers are created at random positions, so that n remains 



This model was first developed by Caroli and Nozieres |3fl| . Al- 
though their formal results directly apply here, let us point out 
that the original physical interpretation - that the elementary 
mechanically unstable units we re t he interasperity microcontacts 
as a whole — was not correct l-Hl . Let us insist again that the 
relevant instabilities do not occur on the micrometric scale but, 
within the inter- asperity junction, on the nanometric scale. 



constant. We assume, for simplicity, that horizontal dis- 
placements are purely one-dimensional. 

x 
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h 



FIG. 12: (a) Toy model of the adhesive junction. X denotes 
the relative displacement of blob A with respect to B in the 
absence of shear compliance (k — 0), x = X + u the real 
relative coordinate of the shear-deformed A/B contact, (b) 
The repulsive interblob potential U(x). 

The system is sheared by imposing to the upper plate 
the displacement X with respect to the lower one. Let k 
be the shear stiffness of the system formed by two con- 
tacting blobs. Under the action of the pinning poten- 
tial, the center of each blob contacting surface undergoes 
an elastic shear displacement u/2 along Ox, so that the 
horizontal separation x = X + u. The total energy of a 
contact is thus : 



£{X,u) = U{X + u) + -ku 2 (34) 

Assume for the moment that T = 0. Then, for fixed 
X, u sets at the value u*{X) which minimizes £, such 
that : 



U'(X + u*) = -ku* (35) 

Two cases must then be distinguished, depending on 
whether the instantaneous equilibrium thus defined is 
unique (monostable contact) or not (bistable contact). 

Monostable contact: 

If k > ko = max(—U") ~ Uq/o 2 , - with a the range 
of U, Uq its maximum value - the solution of eci. (|33|l is 
unique whatever the value of the reference coordinate X 
(see Figure EH- 




FIG. 13: Graphical solution of the local equilibrium equation 
13511 in the monostable case. 
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Assume that we impose a constant external shear ve- 
locity X = V. As X increases, u nearly follows its instan- 
taneous equilibrium value u*(X), and the contact energy 
evolves as £*{X) = £ (X, u*(X)). Indeed, the excess 
elastic energy is dissipated out of the contact region of 
size ~ a via acoustic radiation, at the characteristic rate 
Tg^ 1 ~ c/a, with c a sound velocity, hence much larger 
than the shear rate 7 ~ V/h ~ V/a. 

The instantaneous pinning force exerted by the contact 
on the upper plate is simply : 



fp = ~ 



d£* 
~dX 



d£ 
dX 



d£ 
du 



du* 
~dX 



(36) 



i.e., taking advantage of eq.lJSSJ 



f p (X) = -U'(X+u*(X)) (37) 
and the net work spent to sweep through the contact : 



dXf p (X)= U(X + u*))\^ = (38) 



The total instantaneous friction force F(V) for a junction 
of area E is the sum of the pinning forces f p over the nE 
randomly positioned contacts of transverse range a, so 
that : 



V 1 poo 

F = — — / dX (-f p (X)) = nE- = 
a J -00 a 



(39) 



As could be expected, for an elastically monostable 
junction, frictional dissipation vanishes when V — > 0. 

Bistable contact : 

In the opposite case where k < ko, in a finite range 
of values of X equation l|35|l has three solutions it? (X) 
(Figure I14|) , the two extreme ones u\ 2 corresponding to 
minima E\ 2 (A) of the energy £ (X, u), while Uq is asso- 
ciated with a maximum. As X is slowly increased from 
—00, the instantaneous u* evolves continuously with £, 
i.e. follows branch u* (Figure Hlb). The £* equilib- 
rium gradually evolves from fully to meta-stable, while 
the barrier A£b(X) separating the two basins of £ de- 
creases until the spinodal limit X = X + where it van- 
ishes. At this bifurcation point, branch (1) terminates, 
u\ becomes unstable, and the only choice for the sys- 
tem is to relax towards the, now single, lower minimum 
£2- The corresponding finite energy difference A£\2 is, 
again, lost via acoustic radiation, i.e. instantaneously on 
the scale of the drive. The pinning force f p exhibits the 
hysteresis cycle common to such bistable systems (Fig- 
ure !14b ). The work for sweeping through the contact (the 
hatched area noted (1) on Figure ITlh) w = A£i2 is now 
finite. 




£(x,u) 





FIG. 14: (a) Local equilibrium construction for a bistable 
contact for various values of the reference coordinate X. (b) 
Evolution of the contact energy £(X, u) for increasing X. The 
extrema conrrespond to the intersections labelled in (a), (c) 
Pinning force cycle for the bistable contact. When driven 
towards increasing X, the system follows branch (1) up to 
Ci, where it jumps down to C2. 



The friction force F on the whole junction corresponds 
to a random distribution of X, i.e. to a uniform popula- 
tion P(X)dX = n'SdX/a of branch (1) of the hysteresis 
cycle (Figure [T5k). so that Fd = nT,w/a. Multistability 
results in a threshold rheology and, as long as V <C c, 
the dynamic friction force is velocity independent 13 . 




(a) 



(b) 



13 The corrections resulting from imperfect adiabatic adaptation 
are easily shown to be of relative order (V/c) 2 / 3 l39l 
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FIG. 15: The thick lines correspond to the populated parts 
of the force cycle for (a) the sliding junction (b) the junction 
at rest after turning off the shearing force. A is the associated 
rigid recoil. 



Let us now assume that, starting fom this steady slid- 
ing regime, we suddenly suppress the external shear force. 
The interblob contacts must react so as to bring the sys- 
tem to global mechanical equilibrium, where the sum of 
the pinning forces vanishes. The only way for the junc- 
tion to realize this situation is via a recoil A of the upper 
plate : each individual reference coordinate Xi decreases 
to [Xi — A), yielding a new distribution P eq (X) such that 
the two shaded areas on Figure [T3b be equal. P eq has a 
discontinuity at the Maxwell plateau of the force cycle. 
If, starting from this state, we now gradually increase the 
external shear force F, P eq rigidly shifts to the right by 
6X(F). As long as SX < A , no irreversible "jump" oc- 
curs, the displacement is therefore fully reversible : the 
junction is in an elastic regime. This terminates when 
the discontinuity of the shifted distribution reaches the 
spinodal limit X + , i.e. when SX = A, i.e. F = Fd- In 
this stop and go scenario, the static threshold force is 
equal to the dynamic one. 

However, due to multistability, global equilibrium can 
be realized by a huge number of configurations of the 
interblob contacts. Let us assume that we create the un- 
sheared junction by bringing the two plates into contact. 
Depending on the details of this "mechanical quench" , 
each newly formed interblob contact may correspond to 
a state on either of the two branches of the force cycle, 
the only condition to be satisfied being that : 



dX F(X).f pl {X) = 



(40) 



i=l,2 



where i labels the branches of the force cycle. The elastic 
regime ends when the edge of the shifted P\ reaches X + . 
This point, where irreversible sliding starts, is the static 
threshold. Contrary to dynamic friction, it is not an 
intrinsic property of the junction but depends on its past 
history. If, for example, Pi{X + ) ^ 0, dissipation starts 
at vanishingly small shear. This has been illustrated 
in ref. [41| where it was shown that the dispersion in 
measured values of /i s for a PMMA MCI is strongly 
reduced by preparing the initial state in a controlled 
way, i.e. by sliding before repose. 

Rate effect at finite temperature : 

Consider a bistable contact, swept along branch 1 of 
the force cycle, i.e. sitting at the left minimum of the 
£(u) potential (Figure PHI As soon as T ^ 0, thermal 
noise is able to activate jumps over the barrier A£&(X), 
i.e. to provoke premature jumps onto branch 2 before the 
spinodal limit X + is reached. For fixed X, the jumping 
rate is the Kramers one : 



t(X) 



uj a exp 



A£ b {X) 



k R T 



(41) 



with Lu a an attempt frequency, here typically ~ c/a. As 
X approaches X+, A£i, decreases smoothly to zero, and 
t(X) decreases exponentially up to a cut-off fixed by vis- 
cous (acoustic) dissipation. In practice, only close below 
the spinodal limit does the activated jumping rate be- 
come non negligible. The steady sliding distribution P is 
now controlled by the competition between : 

- the advection imposed by the external drive of X at 
velocity V, 

- the thermally induced jumps from branch 1 to branch 

2. 

This is expressed by the evolution equation 14 : 



dP _ dP_ P 
~dt ~ ~ dX~ t{X) 

the steady solution P st of which reads 



(42) 



P s t{X) — Const exp 



x 



dX' 



(43) 



Due to the exponential variation of r, P st is quasi- 
complctcly depleted between X + and a cut-off X c be- 
low which activation plays a negligible role. The friction 
force is thus smaller than its T = limit. Moreover, the 
larger V, the less time the contact spends in the vicin- 
ity of a given X, the less the probability for it to jump 
prematurely. Hence, X c increases with the driving veloc- 
ity, and the friction force increases accordingly : thermal 
noise results in a velocity- strenghtening rheology. 

This rate dependence can be evaluated analytically as 
long as the the barrier height evolution for X ~ X c can 
be approximated simply. For not too large V, where 
variations are linear, it is found (see ref. 0], Chapter 11) 
that : 



F(V) = F(V ) 



1 + a th In ( — 



(44) 



with Vo is a reference velocity in the above-mentioned 
range of validity, and : 



k B T 

VV act 



(45) 



and a ~ nw/a is on the order of the T = frictional 
stress. At larger velocities, when X c comes very close 



14 This expression assumes implicitly that the drive is slow enough 
for the instantaneous jumping rate to be that for the non- 
advected system 
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below X c , A£b ~ (X c — X) 3 / 2 , and the rate dependence 
of F commutes from log-linear to (In I/) 2 / 3 [33] - a be- 
havior which has been observed by Sills and Overney p3 | 
in a nanoscale friction experiment on glassy polystyrene 
performed with an atomic force microscope. 

b. From the toy model to the real junction : 
When noticing that the functional form predicted by 
the toy model (eq.l^H) for the intcrfacial rheology does 
fit the empirical expression (eq.J2H|)) deduced from ex- 
periments, one is led to go one step further. That is, let 
us take for a moment the toy model at face value and 
compare experimental results for coefficient a fea. (|30l) ) 
with ath fea. lj4"5|0 . As mentioned above, Nakatani 
found that, for granite, a increases quasi-linearly with 
temperature, in agreement with equation i|45() - a re- 
sult confirmed for PMMA in the more restricted velocity 
range investigated in [23 . These authors have then been 
able to deduce values of the volume v ac t- In both cases, 
they find it to be on the order of a few nm 3 . 

This naturally leads to conclude that, inspite of its 
crudeness, the toy model does capture the main physi- 
cal mechanisms responsible for the frictional rheology of 
MCI, namely : 

(i) When sheared, the highly confined junctions be- 
tween asperities are the seat of mechanical instabilities, 
each of which primarily affects a small region containing 
a few atomic (or, for polymer glasses, monomer) units. 

(ii) The strengthening, logarithmic, rate effect can be 
assigned to premature "flips" of these clusters induced 
by thermal noise. 

Feature (%)) is fully consistent with our description of 
the junction as a quasi-2Z? disordered medium, solidi- 
fied under the high confinement conditions imposed by 
the bulks of the contacting asperities, but mechanically 
weaker than these bulks - probably due to a higher frac- 
tion of quenched free volume. Indeed, a number of nu- 
merical studies, pioneered by Argon et al (see also 
Falk and Langer |44| and Malandro and Lacks °f 
sheared glasses, both polymeric and molecular, at tem- 
peratures <C T g , have shown that elastic dissipation in 
these systems occurs essentially via sudden collective re- 
arrangements of clusters of typically nanometric volume. 
These clusters, termed by Falk and Langer shear trans- 
formation zones (STZ), are randomly located 15 and their 
average density is constant at constant shear rate 7 16 : 
on average, once a STZ has flipped into a lower energy 
local equilibrium state, another one appears at an uncor- 
rected position. That is, we may interpret the interblob 
bistable contacts of the toy model as a sketchy repre- 
sentation of STZ's, the volume v ac t providing a rough 
evaluation of the average zone volume. 



15 In the absence of shear localization 

16 To which extent and how this density depends on 7 is an impor- 
tant though still unsettled question. 



Feature (ii) opens onto a much more difficult, still 
essentially unsolved question. In the toy model, the 
bistable units are completely decoupled, and their contri- 
butions to the friction stress are simply additive. The real 
junction is, of course, dense everywhere, so that each STZ 
must be understood as embedded in a quasi-2Z? elastic 
medium, and coupled to the adjacent deformable asper- 
ity bulks which, though stiffer, do transmit stresses. The 
flip of a STZ is therefore akin to the tranformation of 
an Eshelby inclusion ^(|: as already emphasized in [4^ . 
it gives rise to a multipolar force signal which deforms 
the surrounding medium [4j|, hence resulting in loading 
steps on the other STZ's, which were "on their way to- 
wards flipping" . The randomness of these signals in time 
and space is the source of the so-called dynamic noise, 
which acts in parallel with the thermal one. Elastic cou- 
pling being long ranged, dynamic noise is likely to trigger 
cascades (avalanches) of correlated flips, which are not 
accounted for in mean field approximations - hence the 
difficulty of this class of problems . 

Moreover, clearly, the larger the imposed strain rate 7, 
the more frequent the flips hence, roughly speaking, the 
larger the effective strength of the dynamic noise - which 
must vanish in the vanishing 7 limit. So, if these effects 
can be modelled in terms of an effective temperature, this 
must be a growing function of 7. Up to now, no theory 
of this effect is available, and it is therefore not taken 
into account by the recent phenomenological theories of 
plasticity of amorphous media : the STZ-based one |44| 
and the soft glass rheology of Sollich et al. |4q|. Qualita- 
tively, one expects such rheologies to exhibit a crossover 
between two limiting regimes : a thermal noise-controlled 
one at low 7 and/or high T, a dynamic noise-controlled 
one in the opposite conditions. 

Our above analysis of friction leads us to conclude that 
the rheology of MCI junctions is controlled , at and above 
room temperature, by thermal noise, in the investigated 
range of 7 — V/h i.e. typically 7 < 10 5 . 




iq, 



FIG. 16: Schematic representation of Persson's model Q. 
The junction is constituted of a dense set of elastically cou- 
pled nanoblocks driven by a rigid plane through shear springs. 
When the stress on a block a < a a , it is stuck to the lower 
plate. When a reaches a a it gets depinned and performs 
damped motion. When its sliding velocity in vanishes it re- 
sticks. 

Persson has studied numerically a model which 
represents the junction, as sketched on Figure l(To]l . as a 
dense set of elastically coupled pinned blocks driven by 
a rigid plate. He found (see Q, chap. 11) that dynamic 
noise effects were negligibly small, thus justifying a 
thermal noise-controlled mean field rheology (equation 
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(|44(l ) . Note however that the smallness of dynamic noise 
effects in his model is likely to result from the infinite 
stiffness of the confining plates, which should lead to 
exponential screening of elastic couplings. 

Our analysis of MCI rheology can therefore be summed 
up by the artist's view of Figure (|17[) . The "threshold- 
plus-logarithmic" behavior fitted by the toy model results 
should therefore be understood as corresponding to an 
intermediate regime, with limits on both the high and 
low velocity sides. 



thermal noise regime 



dynamic noise regime 




V exp al . range 



FIG. 17: Schematic interfacial stress vs velocity curve (see 
text). The Newtonian Eyring regime corresponds to < V < 
V. 



In view of the current high activity in the field of soft 
glass rheology, one may be optimistic about the emer- 
gence of theoretical predictions about the position of the 
high-y crossover towards the non thermal regime. How- 
ever, it must be pointed out that MCI friction is cer- 
tainly not a good tool for experimental investigation of 
this question since, at the necessary sliding velocities, in 
the mm/sec range or larger, one would have to cope with 
the tricky problems related with e.g. frictional heating. 

The opposite low-V limit is essentially academic. Let 
us come back to the toy model : as the stress decreases, 
the cut-off X c of the contact distribution approaches form 
above the position corresponding to the Maxwell plateau, 
and the rate of back jumps from branch 2 to 1 becomes 
non negligible and comparable with the 1 — » 2 one. At 
the same time, the interbranch barrier A£b increases to 
~ cry v ad, making jumps extremely rare. This results in 
the classical Eyring (sinh) behavior for a s (V), ending in 
an extremely fast quasi-linear drop to zero, with slope 
rj/h, where rj is the viscosity of the glassy junction. The 
crossover velocity V* is in practice much to small to be 
observable. 



D. Sliding dynamics of a MCI : 

1. The Rice-Ruina friction law : 

We can now sum up the above phenomenological anal- 
ysis into a model of MCI dynamical friction which reads, 
with the help of equations lfT7Jl. J22J|,|(23J| and O : 



F((j>,±) = a s0 T, r0 



1 + m In 1 



1 + a In 



= 1 xcj) 
D 



(46) 



(47) 



We saw that both m and a are typically of order 10~ 2 , so 
that it appears reasonable to neglect the ma In 2 term in 
expression <|4t)|) for the friction force. Moreover, we saw 
that the short time cutoff r associated with the early 
stage of area creep becomes relevant only for velocities ~ 
Do/t, on the order of a few hundred /zm/sec at least. For 
the moderately accelerated instationary dynamics which 
we will analyze below, in order of magnitude <f> ~ Dq jx so 
that, in the low velocity regime considered here, (j)/r <C 1. 
Equation (|46|l then reduces to the following expression for 
the dynamic friction coefficient : 



[i d (4>,x) = fJ. d (V ) + -Bin 



D 



A In 



(48) 



where the reference velocity Vq can be chosen at will in 
the above-mentioned velocity range. The constants B 
and A are defined in equations ((19(1 and 1)30(1 . 

Equations |^7| ) and |^<S| ) constitute the Rice-Ruina 
(RR) model, proposed by these authors in 1983 on the 
basis of Dieterich 's experiments. 



2. The RR dynamics of a driven block : 

The question then immediately arises of analyzing and 
testing experimentally the predictions of the model in 
terms of sliding dynamics. First of all, does it correctly 
describe the stick-slip (SS) oscillations of a driven block, 
and their disappearance at high enough stiffness? 

Consider the system depicted on Figure ^ Let us as- 
sume for the moment that the sliding velocity is uniform 
along the interface (homogeneous sliding), i.e. that the 
block has a single degree of freedom, the position x(t) of 
e.g. its center of mass. Its equation of motion reads : 



Mx = -K(x - x (t)) - Wfj. d {^,±) 



(49) 



where fid and 4> are specified by equations l|48|) and l|47|l , 
and (x — xo(t)) is the instantaneous elongation of the 
driving spring. 

At any pulling velocity V there always exists a steady 
sliding solution, namely : 



x = V x(t)-x Q (t) = j^n d Qv) (50) 

In view of the non-linearities of the friction law, one 
must wonder about its dynamic stability, i.e. perform 
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a standard linear stability analysis : setting x = V + 
5x(t), <fi = D /V + 5<f>(t), one linearizes the dynamical 
equations in Sx, 5<f>. Thanks to the time invariance of 
the basic state, the solutions for them are of the form 
Const exp(if2i). Steady sliding is stable (resp. unstable) 
when Trofi > (resp. < 0). 

This calculation is performed in Appendix C. It shows 
that, for given values of M, W and V, steady motion is 
stable for K > K c , with the critical stiffness given by: 



K c ^ = (M0 ~ Mi) 



MV 



WDoUi 



(51) 



At the bifurcation point (A" = K c ) Q, is pure real and has 
the value : 



V / /V - Mi 

Do V Mi 



(52) 



That is, the corresponding bifurcation is of the Hopf 
type, signalling that for K < K c motion should become 
oscillatory. 

In equations /i^ = d^dl d '(In 0), = 

d fj, d/d (In x)V , both derivatives being evaluated at <f> = 
D /V, x — V . The RR expression results in M?> = B, 
fi± = A. In the experiments aiming at characterizing 
the SS bifurcation, the block was sliding under its own 
weight : W = Mg. Then with D ~ l^m/sec, A ~ 1(T 2 , 
the inertial correction in the second factor of the r.h.s. 
of equation (|51|l is of order l0^ 5 (V^ m / s ) 2 . So, inertia 
is negligible for V < 100/im/sec. Then the RR model 
predicts that the critical stiffness 



K C = (B- A) 



W 
A) 



(53) 



should be velocity-independent, while the critical pulsa- 
tion 



v 

A) 



B- A 



A 



(54) 



The SS bifurcation has been studied on symmetric 
MCI involving paper 0, glassy PMMA and PS j2t|. In 
these experiments, K was kept constant. Then, when 
M was increased at constant V, the initially steady mo- 
tion was observed to bifurcate at a value (K/M) c beyond 
which it develops oscillations whose amplitude grows con- 
tinuously as K /M decreases, until true sticking phases 
appear (Figure Tracking the bifurcation at various 
l/'s leads to the stability diagram in the plane (K/M, V) 
of the control parameters, a typical example of which is 
displayed on Figure It is seen that (K/W) c is not 
strictly constant, but decreases slowly with increasing V 
at a rate on the order of 10 per cent per decade. On 
the other hand, V/Q c was measured to be constant up to 



experimental accuracy 17 0- So, the predictions of the 
RR model appear qualitatively satisfactory. 
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FIG. 18: Spring elongation (see inset) versus time for a pa- 
per/paper MCI pulled at the fixed velocity V = 1/im/sec. As 
K/W decreases, sliding motion bifurcates continuously from 
steady (upper curve) to periodic stick-slip of growing ampli- 
tude. 
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FIG. 19: Dynamic stability diagram of a PMMA/PMMA 
MCI in the low velocity range. Steady sliding is stable above 
the bifurcation curve (K/W) C (V). 

The parameters A, B, Dq of the model are then eval- 
uated as follows : measurements of dn d /d(\n V) yield 
(B — A), Do is obtained by equation 1)53(1 from the value 
of (K/W) c at some reference velocity - e.g. l/im/sec 
- then A is determined with the help of equation 154J . 
For example, for a PMMA MCI, one thus obtains [5(j : 
A = 1.2 1(T 2 ± 21CT 3 ,5 = 2.3 lCT 2 ± 2KT 3 , D = 
0.4±0.04/im. Such parameter values are fully compatible 
with those deduced from static aging and velocity jump 
experiments. 

In order to test more thoroughly the validity of the 
phenomenological model, we have analyzed whether it 



17 The poor accuracy on fi c was due in particular to the difficulty 
of extrapolating measurements necessarily performed at finite 
amplitude to the continuous bif urcation while, in this case, non 
linearities develop very fast l4ii . 



19 



correctly predicts a few other salient features of MCI 
block dynamics, among which : 

(i) Non-linear development of SS-like oscilla- 
tions close to the bifurcation: it must be mentioned 
first that the RR model as such (constant \x$ = B and 
U± = A) results in a highly singular behavior for K < K c 
|5 1| - namely unbounded growth of the velocity at finite 
time. This unphysical behavior, in contradiction with 
the observation of a continuous bifurcation, can be 
traced back to the simplifying assumption that fi^ and 
Hx are mere constants, which is also responsible for the 
fact that the model fails to accounts for the non zero 
slope of the K C (V) bifurcation curve. Baumberger et al 
|49j relaxed this assumption by assigning the slope of 
K C (V) to a (In 2 <j>) correction to expression (|48|l for the 
friction coefficient. A standard perturbation expansion 
then results 18 . This extension of the RR constitutive 
law accounted very satisfactorily for the growth of the 
oscillation amplitude and of the frequency shift with 
{K c - K). 

(ii) Creep-like sliding motion following cessation 
of the drive: As can be seen on Figure[S]once the drive 
at velocity V is stopped, the block continues to slide 
while slowing down until it stops at a finite force level. 
Baumberger and Gauthier analyzed the dependence 
of such relaxation curves on the control parameters 
K/W, V. They showed that the RR model predicts 
that the total distance slipped before arrest should be 
^-independent, while experiments show a 100% increase 
on one decade of V. Again, this discrepancy is cured by 
including into equation (|48|) the same In 2 cj> correction 
as for non linear effects close to the bifurcation. An 
example of the resulting fit of relaxation data is shown 
on Figure [501 



FIG. 20: Log-Log plot of reduced slip distance x/Dq vs 
reduced time Vt/Do during creep-like relaxation of a pa- 
per/paper MCI. The system was sliding at V = 1/im/sec 
(open circles) and 10/im/sec (full circles) until the drive was 
stopped at t = 0. The lines are the fits obtained from the RR 
model with a common set of parameters. 

These two examples indicate the limits to be assigned 
to the rate and state model in its simplest version. While, 
clearly, the RR law does capture the essential features of 
MCI friction - namely geometric aging and dissipation 
governed by thermally-assisted mechanical instabilities 
on the nanometer scale - it is insufficient to account for 
finer dynamical features involving finite departures from 
stationary sliding. As described above, phenomenological 
corrections to the expression of /i^ permit to extend the 
validity of the model. However, such extensions must be 
regarded with proper caution : clearly, the non-zero slope 
of the K C (V) bifurcation curve indicates that an exten- 
sion of expression is needed. But this leaves open the 
question of the respective weights of the various possible 
corrections which we have mentioned when deriving ex- 
pression H48|) . Not only is a In 2 cj> term possible, but also 
corrections of the form (lni) 2 and (In <f>) (In x) as well 
as higher order ones. Moreover, corrections associated 
with the short age time cut-off t might become relevant. 
This would mean introducing many more parameters, in 
principle to be fitted from strongly non-linear dynami- 
cal features. As shown by the detailed analysis of the 
resp onse to to large amplitude normal force oscillations 
[53 sucn fits become of more doubtful value as the 
dynamical complexity increases. This points towards the 
fact that formal extensions, or regularizations, of such a 
constitutive law, though they may appear useful, are in 
general difficult to legitimate in detail on a physical ba- 
sis. So they must be considered with a critical eye when 
using them for predictive purposes. 



3. Limitations of block models for extended MCI : 




1 

Vt/D„ 



18 The expansion parameter is found to be e ~ j4/M2 with [12 = 
[d 2 fid/ 9 (In x) 2 ], so that for the RR model (fj,2 = 0) the expansion 
explodes. For PMMA, fi2 is measured to be ~ 10 -3 , hence e ~ 10 
: non linearities develop fast in the SS regime. 



Seismic faults may be represented as spatially extended 
MCI under both normal and tangential loading. Study- 
ing their dynamics then means solving for the elastody- 
namic motion of two "semi-infinite" deformable media in 
frictional contact. It was the need for a realistic con- 
stitutive law for fault friction which motivated the work 
which resulted in the formulation of the RR model. 

Once such a law is available, one is naturally led to 
plugging it into a continuum mechanics description. On 
the other hand, note that, since seismic events amount 
to the nucleation and propagation of interfacial cracks, 
they involve variations of deformation fields down to the 
very small scales relevant close to crack tips. 

This leads us to an important remark resulting from 
the above analysis of the physics involved in the RR law. 
We saw that this law emerges from two averaging pro- 
cesses : 
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(i) the rheology described by the interfacial strength 
a s results from averaging over dissipative events involv- 
ing nanometric regions within junctions; its expression is 
thus valid only on scales much above the nanometer one. 

(ii) the creep growth of the real area of contact is trans- 
lated into an geometric age which is an average over a 
large number of microcontacts. It thus only makes sense 
on a scale much larger than intercontact distances, i.e. 
typically millimetric. 

As already mentioned, this makes the much discussed 
question of the ill-posedness [H^ of the continuous limit 
a formal one. In view of the above remarks, which lead 
to the existence of a physically based small space scale 
cut-off, the most physical recipee for curing mathemati- 
cal ill-posedness seems to be the one proposed by Simoes 
and Martins [5{| : they wash out the UV singularity by 
replacing the local friction relation between interfacial 
normal and tangential stresses (<Jt(x) = — [nd&njx) by a 
non local one (<Jt{x) = —J dx'w(x — x')[/idCn]x'), with 
w a spreading function of finite width. 

In practice, the dynamic complexity due to the non- 
linearities of MCI friction [^| can only be studied nu- 
merically. This is implemented by discretizing the sys- 
tem into elastic blocks of size L x L x L 19 . Rice[57j has 
pointed out that the existence of the stick-slip instability 
imposes a natural upper limit on this block size. Indeed, 
such a discretization implies that the degrees of freedom 
of a block reduce to those of, say, its center of mass, 
internal ones being irrelevant. A block of size L has a 
tangential stiffness Kr, ~ EL, with E an elastic modulus 
of the bulk. It bears the normal load Wl = <J n L 2 , with 
a n the average normal far field stress. It is stable against 
SS if K L D /W L >{B-A), i.e. L < L R with : 



L, 



E 1 

cr „ B — A 



D 



(55) 



Assume that a block has a size Lr < L < 2Lr. Cut it 
into 8 subblocks of size L/2. Each interfacial subblock is 
stable vis-a-vis SS, while the block is not. So the whole 
dynamic complexity on scale L is associated with the 
motion of the relative subblock position, which cannot 
therefore be considered irrelevant. 

The maximum size of a discretization block is 
therefore Lr. From equation (|55|l . with a n /E ~ 10~ 3 , 
B — A ~ 10 -2 , Dq ~ 1-10/im, Lr lies in the meter range. 



19 Choosing blocks much longer in the transverse than in the lon- 
gitudinal direction (i.e., in Burridge-Knopoff-like models, com- 
pressive stiffnesses much larger than shear ones) may lead to dy- 
namical artefacts : for the so-called small events which involve 
only a few blocks, the associated elastic fields only affect depths, 
on the order of their lateral extent, much smaller than the block 
height. It is then illegitimate to neglect internal block degrees of 
freedom. 



III. JUNCTION RHEOLOGY : STRUCTURAL 
AGING/REJUVENATION EFFECTS 

It is clear, at this stage, that the rough-on-rough MCI 
configuration only provides quite an indirect access to 
the analysis of junction rheology. Indeed, information 
about the interfacial strength a s can be extracted from 
friction data only after " deconvoluting" it from geomet- 
ric aging effects, necessarily at the expense of precision. 
This might lead to overlooking finer physical features of 
dissipation within sheared junctions. 

That such might indeed be the case emerged from the 
experimental study by Bureau et al °f the response of 
a rough/rough PMMA MCI to an oscillating shear force : 



F(t) =F + f cos uit 



(56) 



biased about a value Fq below the static threshold. 
For very small /, the MCI responds elastically. As 
/ is increased, the system enters a regime in which 
the oscillatory response is superimposed on a slow, self- 
decelerating, gross sliding motion which corresponds to 
a saturating, finite displacement. This regime prevails in 
a narrow amplitude range where F ma x = Fa + / ^ F s , 
where F s is a static threshold. It was analyzed, on the 
basis of the RR model, as resulting essentially from the 
geometric age dynamics. During each force oscillation, 
the system alternates between sliding, hence rejuvenat- 
ing, when F(t) is close below its maximum, and slowing 
down strongly, thus aging during the rest of the period. 
For / levels such that aging barely wins, gross motion de- 
celerates and the system finally jams. Above a threshold 
/>, rejuvenation wins, resulting in indefinite accelerated 
sliding. 

The RR model accounts satisfactorily for the long 
time dynamics in the jamming regime. However, it ex- 
hibits, for small slid distances on the order of a few 100 
nm, a clear discrepancy with experimental results which 
strongly suggests that another rejuvenation mechanism, 
distinct from the geometric one, hence ignored by the RR 
description, might be at work within the junctions. 



A. Accessing junction rheology directly : suitable 
configurations 

This confirms the interest of studying junction rheol- 
ogy directly, with the help of interfacial configurations 
which either are free of geometric aging or permit to cir- 
cumvent this effect. This can be achieved in three differ- 
ent ways. 



1. Rough-on-flat multicontact interfaces 

So far, we have restricted the definition of MCI to inter- 
faces of macroscopic extent between two rough surfaces. 
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It is then tempting, in order to get direct information 
about junction rheology, to try and realize interfaces a la 
Greenwood which would not undergo the microcontact 
birth and death process responsible, for rough-on-rough 
systems, for the non trivial dependence of geometric age 
on the sliding dynamical history. 

Such a configuration has been realized by Bureau et al 
[58| who studied friction between a rough PMMA slider 
(roughness ~ /xm) and a flat and smooth plate made of 
float glass (roughness ~ nm) 20 . 

The interface is of the GW type but, due to transla- 
tional invariance of the track, the integrity of the con- 
tact population is now preserved upon sliding. The load- 
bearing asperities creep as in a regular MCI but, since 
aging is not interrupted by sliding, the geometric age is 
simply the time T elapsed from the creation of the inter- 
face. The real area of contact E r grows logarithmically 
with T so that, in steady motion at velocity V, the fric- 
tion force F = S r (0 = T).a s (V), becomes an increasing 
function of T, as can be seen on Figure |^ One can 
then take advantage of the slowing down of the loga- 
rithmic growth mode by letting the interface "mature" 
by waiting up to large T's (typically 10 hours), then per- 
forming experiments during a comparatively short period 
(e.g. < lhour) during which S r remains quasi-constant, 
so that friction force variations are directly attributable 
to those of the interfacial strength a s . 
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FIG. 21: Logarithmic growth of the steady sliding friction co- 
efficient at velocity V = 50/im/sec for a rough PMM A/flat 
glass MCI vs time T elapsed since the interface was created. 
Adapted from ref. |5g|. 



2. Junctions in the Surface Force Apparatus 

The surface force apparatus (SFA) was first developed 
|Hflj [riol ] to study adhesion. In a SFA, two atomically pla- 
nar mica surfaces, slightly curved so as to realize a cross- 
cylinder geometry (Figure l2~2l . are brought into contact 



Since glass is much harder than PMMA, no ploughing of the 
track by the slider asperities takes place. 



in the presence of a fluid (the lubricant). The interplate 
distance is decreased down to nanometric values on the 
order of a few molecular sizes. In this boundary lubrica- 
tion regime where the fluid is very highly confined, the 
cylinders deform elastically, leading to a Hertz planar cir- 
cular contact, with radius commonly in the 10/xm range. 




confined <" H 
"lubricant" ~10mm 

FIG. 22: The cross-cylinder geometry used in the SFA (upper 
sketch) is equivalent to a boundary lubricated sphere-plane 
contact (lower sketch). 

The SFA was later modified [6l| so as to allow for shear 
motion and measurement of friction forces. Since the area 
is measured optically, the normal stress borne by the con- 
tact is directly accessible, as well as the shear stress asso- 
ciated with frictional dissipation in the junction formed 
by the lubricant. In many instances (see below) such 
boundary lubrication layers behave as weak disordered 
solids, exhibiting an elastic regime at low shear levels. 



3. Extended soft contacts 

We refer here to the contacts which form between a 
highly compliant (soft) slider with a smooth surface and 
a hard flat smooth track. High slider compliance ensures 
that molecular adhesive contact is realized everywhere 
along the interface, even in the presence of submillimetric 
roughness. 

Such contacts can be achieved in either of two config- 
urations : 

— flat- on- flat: up to now, this geometry has been used 
with sliders made of moulded hydrogels (Young moduli 
in the 1-10 kPa range) |62| |63|. The contact lateral 
extent is usually of order centimeters. 

— ball-on-flat: this is the configuration commonly used 
to study friction at elastomer/glass interfaces [|34[ 

|(5(| . Pressing the ball onto the flat under controlled nor- 
mal load results in the formation of a circular adhesive 
Hertz contact ^3|- Since elastomer Young moduli typi- 
cally lie in the 1-10 MPa range, contact radii in such 
experiments are commonly of order millimeters. 
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B. MCI junctions revisited : 

The results described here were obtained by Bureau 
et al [j^ on MCI formed by rough PMMA sliding on 
smooth (float) glass. The surface state of the glass was 
controlled by grafting onto it a single monolayer of short 
silane molecules, which passivate strongly adhesive sites. 

As already mentioned, such interfaces, when "old 
enough" , can be considered to work at constant area of 
contact S r , thus giving direct access to variations of the 
frictional stress 21 . 



1. Structural aging : 

The central result is illustrated on Figure [23 Namely, 
standard stop-and-go experiments (see §II.B.l) reveal the 
presence of a static friction peak, the amplitude of which 
grows logarithmically with the waiting time t w (Figure 
Eg. 
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FIG. 23: Reduced shear force on a rough PMMA/silanized 
flat glass MCI in a stop and go experiment. For t < 
the block slids steadily (V = 10^im/sec). After the drive 
is stopped (t — 0) the slider slows down and stops at a self- 
selected stress level. At t = 30 s loading is resumed at the 
initial velocity. Adapted from [H^|. 

Contrary to what is the case for rough/rough MCI, this 
peak cannot be associated with geometric aging, which 
is not operative here. On the other hand, the very ex- 
istence of a fi s (tw) larger than the dynamic /id means 
that the interface is the seat of an aging-when-waiting 
vs rejuvenating-when-sliding process, which must there- 
fore necessarily take place within the junctions and affect 
their structure. The existence of this structural aging en- 
tails that <j s is not, as we assumed up to now, a mere 
function of the instantaneous sliding velocity x, but does 
itself depend upon the dynamical history of the interface. 



Moreover, one observes that the structural aging rate 
depends noticeably on the value of the tangential stress 
applied during the stop phase (see Figure l2"4"|) : the lower 
its level, the slower aging is (see figure 5 of reference |58| ) . 

Finally, in many cases - and especially for small wait- 
ing stresses - structural aging becomes observable only 
after a finite, strongly sample-dependent, latency time 
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FIG. 24: Logarithmic growth of the static friction coefficient 
fi s with waiting time t w . Aging at rest under : (•) the self- 
selected arrest shear stress; (o) null stress. From reference 




2. Steady sliding dynamic friction : 

Figure |2U shows a typical example of the velocity de- 
pendence observed for [i^- It exhibits a minimum at 
V = V m i n . The velocity- weakening observed for V < 
Vmin provides another evidence that the structural age 
decreases with V, i.e. that sliding induces rejuvenation. 

For V > V m in, l^d increases quasi- logarithmically 
(fJ-d/fJ'min ~ (1 + ochxiy/Vmin)). That is, we recover 
here the expression of a s (equation Q29(l 1 deduced from 
the analysis of experiments on rough/rough MCI, which 
were used to formulate the RR model. Moreover, the 
reduced slope a ~ 4.1CP 2 has a magnitude compara- 
ble with the value (a ~ 5.I0 -2 ) found for rough/rough 
PMMA junctions. This enables us to conclude that al- 
ready quite close above V m i n , the rejuvenation effect is 
quasi-saturated. 



3. Discussion 

This phenomenology, the fine details of which were 
masked, for rough/rough sytems, by the larger effects of 
geometric aging 22 , appears qualitatively consistent with 



21 As the absolute value of E r is not measurable accurately, the 
absolute value of a a cannot be accessed in this configuration. 



As can be seen on Figure l24l the log-slope of {i s associated with 
structural aging is on the order of a few 10 — 3 , typically from 10 
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FIG. 25: Dynamic friction coefficient fid vs velocity V. From 
reference l5al. 



our description of junctions as confined weak glassy me- 
dia. Indeed, it is closely akin to the behavior of sheared 
bulk soft glassy materials, such as colloidal glasses [^J 
[o8| and pastes jg^. The rheology of these systems has 
also been modelized [ijj [zij as resulting from the 
interplay between : 

— aging at rest, i.e. strengthening via slow relaxation 
down the highly multistable energy landscape character- 
istic of quenched disordered, jammed, media. 

— rejuvenation by motion, which can be seen as reshuf- 
fling the populations of the local minima. 

In the sliding (i.e. plastic) regime, the rate of the lo- 
calized dissipative events invoked in §11. C. 4 should be 
increased by the decrease of the energy barriers due to 
stress induced biasing. The "young", shallow energy 
states which were depopulated as aging pushed the sys- 
tem into "older" , deeper ones, thus can get repopulated : 
the sheared soft glass rejuvenates. 

In the same perspective, one may tentatively interpret 
the above-mentioned dependence of the aging rate on 
waiting stress level as follows. In a stop-and-go exper- 
iment, the "stop" acts as a mechanical quench of the 
strongly rejuvenated formerly sliding state. Aging then 
restarts in a landscape with barriers whose height is de- 
creased by stress-induced biasing, hence accelerates when 
the stress level increases 23 . 



to 4 times smaller than the aging slopes for rough/rough MCI. 
This leads one to conclude that, in the latter case, the growth 
of fj, s , though dominated by the geometric effect, does contain a 
small contribution due to s tructural aging. 

Viasnoff and Lequeux |72| have found that applying to a col- 
loidal glass, aging under stress-free conditions, an oscillating 
shear stress of finite duration does result in a strong pertur- 
bation of the aging process. However, it is important to note 
that, in their experiment, the stress amplitude is larger than the 
yield stress of the material, which certainly results in some non- 
stationary flow. This contrasts with the situation of ref. |5S|. 
where the applied waiting stress does not provoke any sliding. 



C. Boundary Lubrication Junctions 

Over the past twenty years a considerable body of re- 
sults has come out of friction experiments performed in 
the SFA, in the boundary lubrication (BL) regime. It is 
much too vast to be extensively reviewed here. We limit 
ourselves to those salient features which we consider to 
be most relevant in the perspective of this article. 



1. Confinement-induced solidification 

When the "lubricant" - a liquid in the bulk phase at 
the temperature of the experiment - is compressed be- 
tween the two mica surfaces, it is gradually squeezed 
out until its thickness is reduced to a few molecular 
sizes. Then, for materials composed of small, spherical 
or short-chain, molecules (e.g. OMCTS - an inert sil- 
icone - squalane, linear alcanes) one observes [6l| that 
the normal force- us-thickness Fn(D) curve develops os- 
cillations whose amplitude grows rapidly as D decreases 
(see Figure I2rj|) . This behavior can be ascribed to the 
layering of the confined fluid, each successive minimum 
being reached by squeezing out one more full layer. 
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FIG. 26: Measured oscillatory force between two mica sur- 
faces immersed in the liquid OMCTS. The arrows indicate 
inward or outward jumps from unstable to stable positions. 
Inset : peak-to-peak amplitudes of the oscillations as a func- 
tion of D. From ref. |61| . 

Thin enough (typically 3-4 molecular layers) such 
junctions, when sheared, exhibit an elastic response at 
low shear stress levels, until they reach a static thresh- 
old where they yield and start sliding [H O [zl- This 
demonstrates that they have solidified under the effect of 
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the high confinement 24 . This behavior has been con- 
firmed by numerical studies 76]. 

Klein and Kumacheva jz3 found that, in the case of 
OMCTS (and at variance with previous results reported 
in |z3), the liquid-to-solid transition occurs abruptly 
when decreasing the number of layers from 7 to 6. 
Whether or not such abruptness is a general feature has 
not been documented on different systems. 

Finally, it was shown in |74[ that increasing the tem- 
perature makes the confined solid weaker, i.e. leads to a 
decrease of its yield stress. 

It is important to note that the highest normal stresses 
which can be realized in the SFA are on the order of 10 
MPa. That is, they are at least one order of magnitude 
smaller than the levels met in MCI microcontacts, which 
we have seen to be comparable with the yield stress of 
the confining bulk solids. 



2. Structural aging 

Stop and go experiments show that the static thresh- 
old of solidified BL junctions increases with the waiting 
time t w spent at rest. This strengthening, which reveals 
structural aging, has been studied quantitatively on a va- 
riety of lubricants, including hexadecane 1741 . squalane 
[78| [7^ 1 and a star-shaped polymer melt |80j |. Strength- 
ening also manifests itself through the increase with t w 
of the layer elastic stiffness |8J . 
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For some fluids made of short-chain molecules, e.g. hexadecane 
|74| , the solid-like response only buils up after a finite amount of 
shear-induced sliding - suggesting that in this case initial sliding 
helps ordering. 



FIG. 27: Stiction peak heights vs waiting time t 3 for (a) 
hexadecane (from ref. (b) squalane (from ref. L 

is the normal load, t c the latency time. 

In all cases, strengthening is slow, namely either linear 
on a logarithmic t w -scale [74| [glj or somewhat faster [T^ . 
As illustrated on figure |2"71 finite latency times, ranging 
up to about 10 seconds, previous to aging, have been 
observed on hexadecane [24| and several lubricants in- 
cluding squalane [78j]. This result was later contradicted, 
for squalane, in reference [zi| ~~ a discrepancy which is 
possibly ascribable to the difficulty of defining such a no- 
tion with high precision. Indeed, at short waiting times, 
the height of the stiction peak (the difference between 
the static threshold force and the dynamic one at the 
imposed driving velocity) becomes very small. Since ac- 
curacy is limited by the noise level, it may be difficult to 
ascertain whether one should invoke complete latency or 
a gradual increase of the aging rate. 

Finally, Drummond and Israelachvili |78| have com- 
pared, for squalane and a star-shaped lubricant (PAO), 
aging under zero and a finite shear stress close below 
the dynamic level. Contrary to what was observed on 
PMM A/glass MCI, they find that, for a given t w , the 
stiction peak is systematically higher when waiting un- 
der zero shear stress. This leads the authors to associate 
rejuvenation upon sliding, responsible for the force drop 
producing the stiction peak, with shear-induced molecu- 
lar alignment, whose relaxation is certainly faster at zero 
waiting stress. We will come back to this in the discus- 
sion of §HI.F. 

3. Sliding dynamics 

All BL systems with a finite static threshold exhibit 
the same dynamical behavior in the explored low driving 
velocity domain (ranging from < 10 _2 /zm/sec to a few 
10/im/sec). Namely, as illustrated by figure 1251 : 

— For V < V c i periodic stick-slip (SS) motion. 

— In a finite range V c \ < V < V C 2, the motion exhibits 
intermittency, in the form of randomly spaced stick-slip 
peaks HI]. 

— For V > V C 2, these peaks disappear, sliding becomes 
steady |82j ] and, in this regime, dynamic friction is V- 
weakening. 

In most cases studied, the amplitude of the SS force 
signal is roughly constant up to V C 2 ■ 

It is worth mentioning that a qualitatively similar be- 
havior is observed with gelatin/glass extended sliding 
contacts [83|. In this latter case it is associated with 
the fact that, in the SS regime, sliding is spatially in- 
homogeneous - it occurs via propagating slip pulses (see 



Whether or not such a discontinuous transition exhibits hys- 
teresis upon increasing/decreasing V is not documented, to our 
knowledge. 
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suit in the formation of a thin interfacial layer of solvent, 
hence to hydrodynamic lubrication. For attractive cou- 
ples, friction is strongly influenced by the formation of 
transient adhesive polymer-substrate bonds, which they 
modelized by adapting a description due to Schallamach 
(see below). 

— the dependence of the friction force on normal loads 
in these different situations. 

However, up to now, the only system on which the 
nature of the sliding dynamics has been fully analyzed is 
gelatin on glass |63j, on which we now concentrate. 



FIG. 28: Dynamic friction force vs time for a squalane bound- 
ary lubrication layer driven at increasing velocities under pres- 
sures P = 1.4 MPa (a); 5.4 MPa (b) (adapted from Figure 2 
ofref. 

§111. C. 4 below). In the intermittent regime above V c \, 
the corresponding large force peaks are interspeded with 
smaller ones due to nucleation, propagation and death of 
small slip events within the contact. This remark sug- 
gests the possibility that sliding in BL junctions might 
also exhibit spatial complexity. Such a model has been 
proposed by Persson [3j , who invokes the nucleation and 
propagation of "shear- melted islands". 

In an extensive study of squalane junctions, Gourdon 
and Israelachvili |7S| have, in particular, investigated the 
effect of normal load on the sliding dynamics. They find 
that the dynamics described above prevails at high nor- 
mal stresses (P > 5 MPa). However, at lower normal 
stresses ( < 2 MPa), a different behavior is observed : 
for V > V c i, the amplitude of the erratic stick-slip de- 
creases continuously down to zero. 

Finally, using a SFA with a large sliding range 
(500/im), Drummond and Israelachvili 78] have studied, 
for squalane, the shape of friction force transients follow- 
ing velocity jumps. These transients are long-lived, and 
characterized by a distance comparable with the contact 
diameter. 



D. Extended Soft Contacts : Gelatin/Glass 
Friction 

They are particularly easy to realize with hydrogels, 
such as gelatin, which exhibit both ultra-low shear 
moduli (typically a few kPa) and an extended elastic 
regime. Their high solvent content (typically > 80%) 
endows them with specific frictio nal p roperties resulting 
from their poroelastic character j8J]. For this reason, 
we separate this case from that of non-swollen elastomers. 

Gong and Osada 62] have performed extensive studies 
of: 

- the influence on the dynamic friction level of the 
(attractive vs repulsive) nature of gel/substrate interac- 
tions. They have shown that repulsive interactions re- 



1. Static threshold 

The experiments described here were performed under 
zero normal load, adhesion being strong enough to lead 
to solid friction, i.e. to a finite static threshold. 
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FIG. 29: Schematic representation of a gel block driven at 
velocity V > V c (see text). The self-healing slip pulse propa- 
gates at velocity Vu p . 

When a gelatin block is sheared (Figure |25J), as the 
force increases, the gel/glass contact first remains pinned, 
and the block deforms elastically. At a threshold stress 
do sliding sets in. This occurs via nucleation at the trail- 
ing block edge of an interfacial fracture without mea- 
surable vertical opening, which propagates forward with 
a constant velocity Vu p in the mm-cm/sec range, i.e. 
much smaller than that of transverse sound (typically 
~ 1 m/sec) 26 . By changing the gelatin concentration, 
hence the mesh size £, and the solvent viscosity rj s , it 
appears that Vu p ~ D/S,, with D the collective diffusion 
coefficient [84| , measured from quasi-elastic light scatter- 
ing experiments [85| : D ~ G£ 2 /r) s ~ hsT/ris^, with 
G ~ ksT/^ 3 the shear modulus. 

This indicates that the interfacial pinning, which must 
be severed in order to crack, is provided by bonds with 
a typical lateral spacing d on the order of the mesh size. 

Indeed, the value of the self-selected Vu p can be inter- 
preted as follows. The energy released by bond snapping 
events is emitted at the frequency to ~ Vu p /d, and wave 
vector q ~ g? -1 . The lowest deformation mode of the gel 
block is the so-called Biot mode, associated with network 
diffusion within the solvent ojd = Dq 2 . It resonates with 



Note that, while the shear modulus is controlled by the elasticity 
of the loose polymer network, the undrained bulk modulus is that 
of the solvent, so hydrogels can be considered incompressible. 
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emission, hence becomes most efficient to dissipate the 
fracture energy, when lujj ~ to, from which it emerges 
that d ~ £ 2T - 
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FIG. 30: Aging-induced growth of the static friction peak 
(curves 3,4,5) for a 5% gelatin gel/glass interface. If the 
drive is resumed before the interface resticks (left of vertical 
arrow) no peak is observed (curves 1,2). From ref. [g^J. 

Now, one also observes, as shown on figure 1301 that 
the threshold Co for crack nucleation is not a constant, 
but increases logarithmically with waiting time at rest t w 
(in the explored range t w > lsec). This indicates that 
adhesive bonding between the substrate and the polymer 
tails and/or loops dangling from the disordered network 
proceeds on two time scales : a rather fast one, which 
leads to finite static pinning with a bond areal density 
~ £ 2 , followed by a slow, thermally activated relaxation 
of the adsorbed polymer configurations leading to the 
slow increase of this density. 



2. Sliding dynamics 

When the gel block is driven at velocity V, two types 
of dynamics are observed. 

(i) Above a critical velocity V c (ranging, depending 
on gel composition, from ~ 25 to ~ 250/iin/sec), slid- 
ing is steady. The sliding stress u s (V) is velocity- 
strengthening, and corresponds to a shear-thinning rhe- 
ology, namely : 

with a = 0.6 ± 0.07. This is reminiscent of the behavior 
of polymer solutions and suggests that, in this slid- 
ing regime, interfacial pinning plays only a minor role, 



It is worth mentioning that Rubinstein et al 18(1 heve recently 
observed that sliding between purely elastic solids also sets in 
via propagation of interfacial cracks, but in the elastic case these 
propagate at velocities in the sonic range. The low values of Vu p 
for gelatin are therefore, in contrast, a signature of poroelasticity. 



dissipation being due to the viscosity of a layer of thick- 
ness the mesh size (on the order of 10 nm) formed by a 
solution of polymer segments attached to the network. If 
such is the case, the shear strain rate 7 ~ V/£, and the 
effective viscosity 



Veff = -4 ~ Vs(jT r ) 01 
1 

with tr a typical relaxation time for a chain length ~ £ 
in solution. Using for tr the Rouse time one gets, for a 
variety of gelatin concentrations and solvent viscosities, 
the collapse of data shown on Figure |^ which lends 
excellent support to the above interpretation. 
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FIG. 31: Reduced effective viscosity ?? e ///?7s versus Weis- 
senberg number We = ^tr for gelatin gels in water/glycerol 
solvents. Filled symbols : gelatin (wt %) 5(circles), 8 
(squares), 10 (triangles), 15 (diamonds) in pure water. Open 
symbols: 5 wt % gelatin in 21 % (circles) and 42 %(squares) 
glycerol in water. 

(ii) For V < V c , periodic stick-slip prevails. Optical 
observation shows that during the "slip phase" sliding 
is not homogeneous, but proceeds by propagation of self- 
healing pulses (Figurc|2nj • Pulse heads are the previously 
described crack tips, behind which the local slip velocity v 
decreases from the standard quasi-diverging field towards 
the driving value V. However, this regular decrease is 
suddenly interrupted when v reaches a value equal to the 
critical value V c which turns out to coincide with the 
upper limit of the stick-slip regime. At v = V c , resticking 
(healing) occurs quasi-discontinuously. 

When V approaches V c from below, the pulse length 
increases, so that the transition to steady sliding takes 
place via increasing time spacing of pulses of quasi- 
constant amplitude, while being non hysteretic. More- 
over, for V ~ V c , an intermittent behavior, already 
described above (§111. C. 3) is observed - the interpulse 
spacing becomes erratic, and they are interspeded with 
smaller force signals associated to short-lived small prop- 
agating slip events. A more detailed characterization of 
this complex dynamical behavior will certainly be of in- 
terest. 
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3. Rate and state interpretation 

Clearly, as for the previously considered cases, static 
strengthening means interfacial aging at rest, which we 
expect to be associated with rejuvenation, i.e. weak- 
ening, upon sliding, as proved by the reproducibility of 
ao(t w ) in stop and go experiments. 

This is consistent with the abruptness of the restick- 
ing process at the trailing edge of the self-healing pulses. 
Indeed, assume that the steady-sliding characteristics is 
V- weakening for V < V m in- Close behind the pulse head, 
the local v, larger than V m i n , lies on the locally stable 
branch of a s (V). As the distance behind the head in- 
creases, v decreases. When it reaches V m i n , the interface 
becomes unstable on a small spatial scale (smaller than 
the optical resolution). In other words, since the driv- 
ing stiffness, provided by the very compliant gel block, 
is extremely low, the interface jumps down to V = at 
quasi-constant stress. This leads to concluding that the 
resticking velocity V c corresponds to a minimum of cr s (V) 
above which, as appeared in the previous paragraph, in- 
terfacial pinning should be negligible. This interpretation 
is also consistent with the coincidence of the resticking 
velocity and the disappearance of stick-slip as soon as 
the asymptotic value V of the local velocity v, needed 
for steady sliding to get established, reaches the stable 
branch of a s (V). 

So, here again, it seems natural to try and define a 
structural age of the gel/glass interface. Quite clearly, 
the corresponding dynamic state variable must somehow 
measure the pinning energy associated with the forma- 
tion of polymer-glass adhesive bonds. 

Already 50 years ago, Schallamach [64| proposed, in 
the context of dry rubber friction, a seminal model of 
such a dynamic state variable, which remains the basis 
of all more recent theories of gel and rubber [s^l HH 
friction. In this model, briefly summarized an discussed 
in Appendix D, the elementary mechanical instability 
(see §11. C. 4) is the snapping of a bonded molecule out of 
its adsorption site under shear loading through an elas- 
tic spring (here the polymer segment). Once this energy 
barrier has been jumped over, under the combined effect 
of loading and thermal activation, the molecule is freely 
advected until it readsorbs via a thermally activated pro- 
cess. Of course, readsorption is thwarted by advection, 
which limits the time available for this process to take 
place. The state variable is thus the density of bonded 
molecules. The friction stress is the product of this V- 
weakening density, and of the ^-strengthening average 
depinning force (the faster advection, the less time for 
escaping above high barriers) . So it exhibits a bell shape 
with a maximum at V — V max . The low velocity regime 
is an Eyring viscous one, with the bond density close 
to its equilibrium saturation level. In the I^-weakening 
high velocity regime the depinning force saturates, cr s (V) 
is then fully controlled by dynamic rejuvenation, hence 
vanishes for V — > oo. 

Models of the Schallamach type call for two important 



remarks : 

(i) The extent of the Eyring-like low V regime (V < V, 
see Figure l55|) crucially depends on the height of the des- 
orption barrier. When, as often assumed - explicitly or 
not - it is on the order of a few fcgT, V is non-ncgligiblc. 
A strong consequence is that, in such cases, no static 
threshold should be observable except if one could load 
at exceedingly large velocities. This points towards the 
interest of systematically associating studies of dynamic 
friction in such systems with stop and go experiments. 

Gelatin/glass interfaces exhibit well defined static force 
peaks and crack tips. We must therefore conclude that, 
in this case, we are dealing with strong bonds - in which 
case V becomes negligibly small (see end of §11. C. 4). 

(ii) If we stick to the Schallamach prediction, in the 
case of gelatin, u s (V) should be uniformly V- weakening, 
since in the explored range V 3> V rnax . How can one 
then explain the existence of the V > V c strengthening 
behavior? 

We meet here with an important shortcoming of the 
Schallamach model. Namely, unpinned molecules are as- 
sumed to glide freely, that is, viscous dissipation in the 
interfacial layer is omitted - although, since its thickness 
is typically nanometric, shear rate levels are high. In the 
case of gelatin we saw that, on the contrary, it is this 
contribution which fully accounts for the observed non- 
newtonian friction at the essentially depinned interface. 
This illustrates the importance of taking into account in 
this kind of problem, not only interfacial pinning, but, as 
well, junction viscosity. When both effects are of com- 
parable importance, this open problem becomes a true 
nanofluidics one, as far as the size of the confined poly- 
mer chains is precisely comparable with the thickness of 
the "channel" , i.e. the junction. 

In summary, once complemented with standard vis- 
cous dissipation in the sliding junction, the Schallamach 
picture, which identifies the state variable with the den- 
sity of adsorbed bonds, provides a sound basis on which 
to modelize gel friction. Clearly, a more elaborate mod- 
elization of the dangling chains, taking into account the 
existence of a multiplicity of possible bonding configura- 
tions |90l ] and a detailed description of chain dynamics 
|88| would be needed to make such models more quanti- 
tative, as well as to touch upon the question of long term 
aging. 



E. Extended Soft Contacts : Dry Elastomer 
Friction 

Due to its importance in a number of applied fields 
(e.g., to name only a few, traction of tires, windshield 
wiper efficiency and durability), this question has been 
a subject of very active investigation in the past six 
decades. We restrict ourselves here to a sketchy summary 
of the commonly accepted concepts used for interpreta- 
tion of interface geometry vs molecular adhesion effects. 
For extensive bibliographies, see [9l| |o6| . 
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We consider here the case of elastomers in extended 
contact with much stiffer solids, which can thus be as- 
sumed non deformable. Bulk viscoelastic dissipation is 
very important in rubbers. Its strong temperature de- 
pendence follows the so-called WLF time-temperature 
superposition principle |9^ |. Grosch [gij], and Ludema 
and Tabor [93, have shown that one must distinguish 
between contributions to the sliding friction stress gov- 
erned respectively by bulk and interfacial dissipation. 



1. Bulk dissipation 

It comes into play as soon as the stiff partner surface 
exhibits non negligible roughness - on a scale larger than 
the average cross-link spacing. Assume for a moment 
that the roughness has a characteristic wavelength A. 
Then, in stationary sliding, a given material point within 
the elastomer feels a stress field modulated at a frequency 
u> ~ V/X, which penetrates down to a depth ~ A below 
the surface. This results in a dissipation governed by the 
loss modulus of the bulk, G"(oj), i.e. obeying a V — T 
superposition principle. One is thus able to account for 
the most salient feature of rubber friction, namely the 
existence of a maximum of the cr s (V) curve, at a veloc- 
ity V m , which increases with the temperature T, above 
which stick-slip motion sets in j94[. 



temperature superposition argument an energy scale ~ 
25 kJ.mol -1 , i.e. ~ 0.3 eV/bond. Whether this energy, 
much larger than the van der Waals ones expected for 
PDMS/silane interactions, is relevant to adhesive bond 
snapping in this system remains an open question. 
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FIG. 32: Temperature-dependent friction stress data for 
PDMS on a silanized silicon surface shifted to room temper- 
ature using a WLF shift factor. Open, gray and black circles 
represent data at 298, 318 and 348 K respectively. From ref. 



2. Interfacial dissipation 

Grosch [94| identified, in his experimental results on 
friction of several elastomers on silicon carbide paper, an 
additional feature - a shoulder on the cr s (V) curve at a 
velocity V m ~ w m A, where u m corresponds to the max- 
imum of the loss modulus and A ~ 6nm. This shoulder 
turned into a broad maximum when sliding on smooth 
"wavy" glass. He attributed this to interfacial dissipation 
- i.e. to the snapping out of advected adhesive bonds as 
modelized by Schallamach (see §111. D. 3). 

It is only recently that progress in surface control en- 
abled Vorvolakos and Chaudhury (VC) [6(| to reconsider 
this question. They used two model systems, consti- 
tuted of cross- linked PDMS sliding on (i) a self- assembled 
silane monolayer supported by a Si wafer (ii) a thin film 
of glassy PS. They were thus able to get rid of surface 
roughness effects, hence of bulk dissipation. 

Silane coverage results in very low energy surfaces, 
so one expects polymer substrate adhesive bonds to be 
rather weak in this case and, hence, Schallamach's model 
to be appropriate. Indeed, VC show that their results 
for u s {V) in this system (see figure l3*2*|) . which extend 
over 5 velocity decades, are well explained qualitatively 
in this frame 28 . Moerover, they deduce from a time- 



PDMS bonding to PS is likely to be stronger than 
that on silane. Indeed, consideration of the VC results 
in the low velocity limit leads to raising the question of 
the possible existence of a finite static friction threshold. 
Systematic investigation of the transient behavior at the 
onset of sliding would certainly be enlightening. 



Finally, it is worth mentioning that non homogeneous 
sliding has been observed [96| |9J at some smooth soft 
rubber/glass interfaces. This occurs via propagating self 
healing pulses. But, contrary to the gelatin case where 
contact is maintained in the sliding regions, these Schal- 
lamach waves are of the mode-I fracture type. That 
is, they consist of regions of interfacial detachment with 
macroscopic opening which re-adhere at the back edge. 
So, in this regime, dissipation is essentially controlled by 
adhesion hysteresis 65] . Though high compliance and 
large viscoelastic losses are thought to promote this slid- 
ing mode, no prediction relating material properties to 
the dynamics, nor even to the occurrence of Schallamach 
waves is available up to now. 



Interfacial viscous dissipation is neglected in the VC analysis. 



As already mentioned, it should come into play, especially in the 
V- weakening regime V > Vm, where it is certainly needed to 
interpret the observed stick-slip behavior. 
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F. A tentative classification : Jammed junction 
plasticity vs adsorption controlled dynamics 

The various frictional behaviors briefly described above 
pertain to junctions which differ in several respects, 
namely : 

• Normal stress level, i.e. confinement pressure P, 
which ranges from zero (gelatin) to about the yield stress 
of the softer bulk material (rough on flat MCI) . 

• Level of adhesive interactions, i.e. strength of the 
pinning sites provided by the confining potential. This 
ranges from very weak (fully silanized substrates) to 
strong (e.g. bare glass). 

• Density of adsorbable units, ranging from dense 
(polymer glasses) to (semi)dilute (hydrogels). 

All these systems share two properties : 

- they exhibit a static threshold 29 

- this threshold ages slowly with time spent at rest. 
The experimental signature of this structural aging is 

the slow (in general quasi-logarithmic) growth of the up- 
per limit of the elastic regime. In experiments where the 
slider is remotely driven at constant velocity V, it results 
in general in a force peak. This unambiguously reveals 
rejuvenation upon sliding, since the sliding velocity v is 
equal to the driving one V both in the steady sliding 
state and at the peak, where the slid distance is minute, 
so that the age is basically given by the waiting time be- 
fore pulling 30 . This means that the state of the sliding 
junction is specified, not only by its instantaneous veloc- 
ity, i.e. by the shear rate 7 it experiences, but also by 
the value of some dynamical state variable. 

On the other hand, different classes of systems ex- 
hibit different steady state rheologies, from logarithmic 
strengthening for highly confined PMM A/glass to power 
law for gelatin. 

While these behaviors call for rate and state descrip- 
tions, it is clear that various classes of systems must 
be distinguished, in terms of the microscopic nature of 
the state variable tp, which is necessarily related with 
the junction structure and, hence, with the pinning 
mechanism giving rise to static friction. The existing 
body of experimental results discussed in §111. B-E is 
still far from sufficient to provide a firm basis for such a 
phenomenology. So, the classification which we propose 
here is only tentative, and primarily intended to help 
structuring further questions to be, in a first stage, 
investigated experimentally. 

Let us start from the very high confinement regime, 
realized in the MCI configuration (§111. B) where, for 



PMMA, the confining pressure P ~ 100 MPa is compa- 
rable with the bulk yield stress. We saw that junctions 
between such a polymer glass and highly silanized glass 
behave as soft glassy nanometer-thick media. 

In other words, we claim that, in such a case, the dis- 
ordered material of the dense junction jams i.e. so- 
lidifies under confinement. It ages as a glass, by relax- 
ing towards deeper local minima of the energy landscape 
(inherent states) and, when it flows, dissipation can be 
attributed to irreversible fli ps o f molecular-sized shear 
transformation zones (STZ) |44j . In this case, frictional 
sliding is nothing but plastic flow of this interfacial glass 
which, being weaker than the bulk, naturally localizes 
shear. 

Now, note that, in the MCI experiments, when the 
glass substrate is well silanized, no direct manifestation 
of the polymeric nature of the junctions is observed. 
However, when the same experiments are performed on 
a poorly silanized substrate |58j , the previously observed 
rather narrow aging peak persists, but is followed by a 
long bumpy stress transient (see Figure 123 ■ Moreover if, 
after reaching the steady state, one performs stop-and-go 
tests with a high stopping shear stress level, the long 
transient is absent when restarting motion. If, on the 
contrary, the stop is performed under zero shear stress, 
the wide bump reappears. This strongly suggests that 
these transients, associated with the presence of stronger 
corrugations of the confining substrate potential, are due 
to gradual stretching and alignment of the constitutive 
molecules, which persists when the shear stress is main- 
tained, and relaxes otherwise. Such rubbing-induced 
alignment has been documented already long ago by 
Pooley and Tabor j9]J on teflon, a material for which it 
is particularly strong. It would be desirable, following 
these authors, to test this interpretation by devising 
a set-up which would permit to change the sliding 
direction by e.g. 90°. 
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Except possibly for some elastomers under low normal stress and 

in contact with highly passivated substrates. 

When established sliding is not steady, but occurs via stick-slip, 

static aging results in the growth of the first peak with stopping 

time. 



FIG. 33: Reduced friction force versus time after a 300s stop 
under zero shear stress. The rough PMMA slider is driven at 
V — 10/im/s on a poorly silanized track. 

This observation naturally raises the question of the 
competition between jamming and adhesive pinning to the 
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substrate. It is reasonable to expect that, the lower the 
confining pressure P, the weaker the jamming, i.e. the 
relative contribution of STZ flips to dissipation. 

In boundary lubrication experiments performed in the 
SFA, pressure levels are smaller that those in MCI junc- 
tions by at least two orders of magnitude. As mentioned 
in §111. C. 2, the amplitude of the transients observed on 
squalane by Drummond and Israelachvili [82j is larger 
when stopping under smaller shear stress. They also 
seem to be associated with quite broad stress bumps. 
This leads one, following the authors, to associate these 
transients with molecular alignment. Whether or not a 
narrow glassy-like aging peak is observable at the highest 
SFA pressures would be worth investigating in detail. 

This remark points toward the need for trying to 
bridge, for given intcrfacial materials, between SFA 
and MCI pressure levels, for example with the help 
of the ball-on-flat configuration. Such comparisons 
are necessary in order to appreciate to which extent 
SFA data can be safely extrapolated to the conditions 
prevailing in macroscopic contacts between hard solids. 

The other limiting situation of purely adhesive pin- 
ning. It is exemplified by gelatin gels in contact with 
glass under zero or very small |63| normal load (§111. D). 
In this case the junction is constituted of a (semi) di- 
lute solution of protein segments in water - no jam- 
ming can be invoked. The pinning responsible for the 
finite static threshold is due to adsorption of the poly- 
mer molecules onto the glass substrate, and one may, 
at least semi-quantitatively, identify the state variable 
ip with the areal density of adsorption bonds. Such a 
structure is able to relax by a combination of thermal 
desorption, chain configuration rearrangements and read- 
sorption which leads, at rest, to the slow increase of ip. 
The number of candidate sites on a segment of course 
depends on the chemical nature of the network-forming 
polymer chains. This dynamical problem, although non 
trivial, might be amenable to a theoretical treatment 
which would extend equilibrium studies such as those of 
Subbotin et al. [9(j along the lines of the Schallamach- 
like steady state model built by Charitat and Joanny [g^ 
for the case where only the end monomer is adsorbable. 

In such a completely non-jammed junction, frictional 
dissipation contains two contributions due respectively 
to : 

— Thermally helped advection-induced depinning. It 
is this mechanism which controls the solid friction behav- 
ior proper for such adsorption- controlled interfaces. 

- Viscous flow of the polymer-water mixture in the 
sheared interfacial layer. 

While the latter is velocity-strengthening, the former 
decreases with decreasing ip, so it is basically velocity- 
weakening. 

Elastomers probably correspond to an intermediate sit- 
uation. The interfacial junctions which they form cer- 
tainly have a dense polymer content, while clearly not 



being solid (jammed), since their bulk itself is not. Junc- 
tion thickness scales as the distance between entangle- 
ments. One of the open questions is that of evaluating 
realistically the viscous contribution to friction of such 
sheared layers. Whether they exhibit, as might be ex- 
pected, slow aging at rest is not documented yet, as far 
as all experiments on elastomer friction have focussed up 
to now on the steady sliding behavior and on its temper- 
ature dependence. 



IV. CONCLUSION 

The route we have chosen to follow in this review has 
led us from macroscopic dynamics down to gradually 
decreasing space scales. Solid friction thus appears as a 
particularly favorable area of material science in which 
the governing mechanisms are "universal" enough to 
allow for a rather general phenomenology down to the 
nanoscopic level. 

A first level of analysis, relevant to the widespread 
case of macroscopic interfaces between rough hard solids, 
identifies two main physical ingredients : 

(i) A shear- induced depinning process, in which me- 
chanical instabilities lead to structural rearrangements of 
nanometric shear transformation zones within the junc- 
tions where shear localizes. It is the associated multista- 
bility of these glassy junctions which is responsible for 
the existence of finite friction thresholds. 

(ii) Geometric aging, i.e. slow creep growth of the area 
of the sparse microcontacts forming such interfaces under 
the high contact pressures associated with this geometry. 
Sliding limits the duration of contact life, making geo- 
metric age a dynamical variable with a memory of the 
sliding history. 

This description provides the physical basis of the Rice- 
Ruina constitutive laws which, when properly extended, 
account for all the main features of the low velocity fric- 
tional dynamics. It also permits to point more precisely 
the limits of such a constitutive law. Among these, an 
important one is concerned with higher velocities — typ- 
ically > a few 100/mi/sec - at which geometric age be- 
comes small, leading to the saturation of the destabiliz- 
ing V-weakening behvior of dynamic friction, and to its 
inversion into V-strengthening. Behaviors in the inter- 
mediate (mm/sec to cm/sec) range are still insufficiently 
documented experimentally to allow for detailed models 
of dissipation in this regime. 

For still faster sliding, self-heating becomes non 
negligible, turning temperature into a new relevant state 
variable which affects junction rheology. An example 
of such an effect is that of the stick-slip dynamics of 
the bowed violin string, which has been modelized by 
Smith and Woodhouse (9^| in terms of the triggering 
by slip-induced heating the of a transition from sticking 
solid to viscous fluid of the rosin rubbed on the bow. 
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Beyond the slow dynamics leading to the notion of ge- 
ometric age, another one emerges, associated with struc- 
tural aging of the junctions, whose effects are masked, 
for rough-on-rough systems, by the larger geometric ag- 
ing ones. Structural aging at rest and its counterpart - 
rejuvenation when sliding - is at work in all types of fric- 
tional junctions, leading to slow growth of static thresh- 
olds with time and to stick-slip-like dynamical instabili- 
ties. It affects the topography of the multistable energy 
landscapes which govern frictional rheology. Further ex- 
perimental information will be needed in order to try and 
specify it in terms of precise state variables, whose nature 
is probably not unique. 

On the basis of comparisons between different classes 
of interfaces, we have suggested a schematic classification 
in terms of two limiting types. 

(i) Jammed junctions, solidified under confinement into 
a soft glass structure. They are relevant at high contact 
pressures such as met with multicontact interfaces. Mul- 
tistability results from their glassy structure itself, and 
in this case frictional dissipation enters the wider class of 
problems constituted of plasticity of amorphous media 
and soft glass rheology. 

(ii) Purely adhesive junctions for which pinning is due 
to adsorption bonds between the junction molecules and 
the substrates (s). They form at contacts under low con- 
finement pressures involving gels and, more generally, 
polymers. 

Many real junctions probably are a compound of these 
two ideal types, their frictional dissipation containing 
contributions of both types whose respective weights 
depend on pressure, temperature, density and on the 
physico-chemical nature of the confining surfaces. 

While, in our opinion, the question of geometric age 
effects in multicontact solid friction is now reasonably 
cleared up, that of juction rheology remains a largely 
open issue worth of further investigation. In order for 
modelization to make significant progress, more detailed 
experimental information is needed, in particular about : 

- Systematic trends associated with varying pressure 
and temperature, for systems with controlled surfaces. 

- The effect of various aging histories, paralleling anal- 
ogous studies in structural and spin glasses. In the case of 
adsorption-controlled junctions, information about aging 
might also be obtained via adhesion tests [99|, provided 
that these be performed on materials with negligible vis- 
coelasticity. 

Numerical simulations of sheared junctions in the 
presence of realistic confining potentials should also be 
of great help. 

Finally, the reader may have been surprised that we 
never used terms such as "nanofriction" or "nanotribol- 
ogy" . This we did purposefully, since we consider them 
rather misleading. Indeed, "nanotribology" is often asso- 
ciated with investigations of the nature of the basic dis- 
sipative processes which, as we saw, occur on the nano- 



metric scale. However, friction, as defined in terms of a 
stress or a friction coefficient, results from the self averag- 
ing of such events within contacts of, at least, mesoscopic 
(micromctric) size. This is precisely why friction laws do 
not make sense on the nanomctric scale. 

We consider that the term nanofriction should be re- 
served for situations where lateral contact sizes are them- 
selves of nanometric order, in which case different effects 
may emerge, s uch as exemplified by aging due to capillary 
condensation |100| . It is also likely that , in such regimes, 
the nature of dissipative processes |101| is different from 
that in mesoscopic contacts. 
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Appendix A : Shear stiffness of a MCI 

Already long ago, a number of studies have tried to 
get independent information about the structure of MCI 
from the load dependence of various interfacial proper- 
ties besi des f riction. These include electrical [Tl| and 
thermal Il02j | co nductances, variations of closure with 
normal load |l03j| , acoustic reflection and inter facial nor- 
mal and transverse mechanical stiffnesses |l04j . In order 
to confront experimental results with the Greenwood- 
Williamson model, each of these properties must be cal- 
culated within this same frame. 

Let us sketch out, as an example, the evaluation of the 
interfacial shear stiffness. For a single Hertz contact of 
radius a between two identical spheres of shear modulus 
G and Poisson ratio v, it was calculated by Mindlin |l05j | 
[13 to be : 



k s — 



AGa 



(57) 



which expresses the fact that the strain energy due to 
shearing is essentially localized, in each medium, within 
a volume of section ~ a 2 and height ~ a. So, the relative 
shear displacement involved in the definition of k s is that 
between any two points on both sides of the contact at 
distances 3> a from it. 

Consider a set of n GW microcontacts between two 
blocks of thickness D, shear modulus G. One may sepa- 
rate this elastic system into an interfacial layer of thick- 
ness ^> a but <C D sandwiched between two bulk regions 
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of thickness ~ D. As far as E# <§; E aj , p , the microcon- 
tacts are elastically independent, and the total interfacial 
stiffness K s is simply the sum of the individual ones, i.e., 
with the notations of § II. C : 



K s = Y—^a (58) 
From equations © , © 



Let us consider a MCI between a linear viscoelastic, 
incompressible and isotropic medium and a rigid one - 
a good approximation for contact between a rough elas- 
tomer and a rough hard solid. We want to describe this 
MCI in the the frame of the simplest version of the GW 
model (exponential distribution of summit heights). 

In a linear viscoelastic material, shear stress and strain 
are related via a retarded elastic shear modulus G(t), by 
the constitutive equation : 



K. 



4G 



W 



w 



{2-v)y/^E s (2 -v){l + v) s 



(59) 



So, the GW model predicts, contrary to immediate in- 
tuition, that the interfacial stiffness of an elastic MCI is 
independent of the bulk material shear modulus. 

Of course, experiments measure the compound re- 
sponse of the interface and the bulk. A more important 
warning comes from the fact that expression (|57|l is valid 
only in the limit of very small ratios of shear to normal 
forces. Indeed, Mindlin has shown that increasing 
shear induces, at the periphery of the contact, a slipping 
annulus along which, would it remain non-slipping, the 
friction sliding threshold would be overcome. This effect 
results in a non-linea r wea kening behavior of k s [l3l | and 
of the MCI response [l06| . 

Consider now a GW interface in the fully plastic limit. 
When a small shear force is applied for the first time, 
it induces some further plastic flow, probably leading to 
some small adaptation of contact radii. Under subse- 
quent unloading-reloading cycles, these adapted contacts 
respond elastically, with an individual stiffness ~ Ga. K s 
is then immediately evaluated to be ~ W/X, where the 
length A is of the order of s/tp, with ip the plasticity index 
defined in section II. C 

The experiments of Berthoud and Baumberger [17] on 
PMMA (a polymer glass)/PMMA and AU4G (an alu- 
minum alloy)/AU4G have confirmed the Amontons-like 
behavior l|59|l of K s . The lengths A for the two systems 
were different, but both on the order of a few microns. 
They were found to be compatible, in order of magnitude, 
with the values predicted from the above model. These 
results confirm the robustness of the three properties of 
MCI, listed in section II. C, which emerge from the GW 
model. 

In view of the difficulty of determining accurately s 
and, even more, R, as well as of the numerous approxi- 
mations involved in the GW model and, even more, of its 
elasto-plastic extension, it would be fallacious to expect 
a more quantitative agreement. 

Amontons-like behaviors have also been measured to 
hold, with various degrees of accuracy depending on the 
specific difficulties inherent to each experimental method, 
on the other interfacial properties listed at the beginning 
of this Appendix. 

Appendix B : The viscoelastic GW model 

Our derivation summarizes the work of Hui et al l.3ll . 



a(t) 



v ' dt' 



(60) 



which can be formally inverted into 



e(t) 



f 

J — < 



dt' J(t - t') 



da(t') 
dt' 



(61) 



J(t), which characterizes the strain response to a unit 
stress step applied at t = 0, increases with time, as a 
result of the progressive relaxation of the constitutive 
polymer molecules. 

Let us first consider the viscoelastic analogue of the 
single Hertz contact (§ II. A. 2), created by applying the 
load w at time t = 0. It is intuitively clear that, since the 
compliance increases, the compression S(t) and the con- 
tact radius a(t) increase monotonously. Lee and Radok 
|107| have shown (see also 01) that in such a case the 
solution of the contact problem is that deduced from the 
Hertz one by simply replacing the shear modulus by its 
retarded counterpart G(t). So, the Hertz load-radius re- 
lation (eq. ©) becomes : 



10 



G(t - 1') 



dt' 



irR 



■a 3 (t>) 



(62) 



with a 2 (t') = R*5(t'). 

We now turn to the corresponding GW problem, con- 
tact under the normal load W being first established at 
t = 0. As time increases, clearly, the interfacial separa- 
tion d(t) will decrease monotonously, the number of mi- 
crocontacts n(t) and their individual area increase, lead- 
ing to the increase of the real area of contact E r (i). 

The three GW equations (eqs.JSJ) - Q) become : 



n(t) — N dz- e 
Jd{t) s 



-z/s 



E r (t) = N / dznR(z - d(t))- e 
Jd(t) s 



(63) 



(64) 
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Wit) = N dz- er z l s / dt' 9{t' - t (z))G(t -t')^-[ -R 1 ' 2 (z ~ d{t')f' 2 



d(t) 



dt' 



(65) 



where 9 is the unit step function, W(t) = W6(t), and 
the function to(z) is the inverse of d(t), i.e. 



In the RR model 



d(t (z)) = z 



(66) 



Equation (|65|l is based upon relation (|62ll and states 
that each contact with summit height z existing at time 
t has contributed since t' = to(z) when it first appeared, 
i.e. for t' > t (z) or equivalently, d(t') < z. 

Setting y — (z — d{t')) / s and performing the y- 
integration in ea. l65|) yields : 



fU = Vdo + B In — + A In — 



(73) 



with Vo some reference velocity, <f>Q = Dq/Vq. However, 
in view of the shortcomings of expression (|73() discussed 
in §11. D, we will leave the functional form of /id unspeci- 
fied at this stage. 

Equations iJHJ, (|72() have the stationary solution : 



W(t) = N (7rRs 3 ) 1/2 f dt'G(t-t')-^(e- d ^/ s ) 
J ~oo at \ I 

(67) 

which can be inverted, with the help of definitions (16011 . 
(ED into : 



,-d(t)/s 



dt' J(t - t') 



,,dW{t') 



On the other hand, from equation l|64(l 

S r (t) = NirRe- d(t)/s 



(68) 



(69) 



so that, finally, since dW/dt 1 — S(t'), the real area of 
contact evolves according to : 



W ( D 
x s t{t) = Vt-a - —\i d I —,V 



V 



(74) 



(75) 



In order to study its (linear) stability, we linearize 
equations iffl). (2) in 5x(t) = x{t) - x st (t);S(f>(t) = 
4>{t) - (j) Q . Then : 



— 5£ H Sx + —Si H — < 

W W V D Q 



(76) 



E r (t) = J—j(t)W 
V s 



(70) 



Appendix C : The driven block : Linear stability 
analysis 

Consider the system defined in §11. D. 2, made of a block 
of mass M, carrying the normal load W , driven at the 
constant velocity V through a spring of stiffness K 31 . 
The equations describing the motion read : 



Mi 



-K(x-Vt-a)-Wp, d (^,x) 



(71) 



V 



i 



-Si = 



(77) 



with/i^ = dfj,d/d(Ln4>), iu,& — dfid/d(lnx)V , both quan- 
tities being evaluated at 4> = Dq/V, i = V. Note that 
we saw that, for slowly sliding MCI, \i± > 0, < 0, 
and /it/y — fi± > 0. The solutions of this system are of the 
form : 



Sx 



(78) 



(72) 



and the frequencies fl of these eigenmodes are the roots 
of : 
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31 We assume that, as is commonly the case, K is much smaller than 
the interfacial stiffness. Finite intcrfacial compliance effects are 



evaluated in 



J 



.M 3 ( m 
W ' 



MV 



V WD 



n 2 



K 
W 



r 



Mi — M<#> 

Do 



in 



K V 

wW 



(79) 



If, for all roots, ImQ, > 0, all (infinitesimal) fluctua- 
tions about the stationary state decay, steady sliding is 
(locally) stable. For very large K the solutions of l|79() 
read : 



(80) 



Since fi± > 0, the system is stable in the large K limit. 
As the stiffness decreases, it becomes unstable if, and 
when, the imaginary part of one at least of the roots 
vanishes. One easily checks on equation Ij79(l that this 
occurs when : 



K 
W 



(m</> - Mi) 1 + 



MV 
WD of n 



(81) 



As expected, instability only occurs when [1$ — \x± > 
0, i.e. (see §11) when the steady sliding fid is velocity- 
weakening, which is indeed the case for our systems in 
the V-range under consideration. 

Note that inertia only comes into play in eq. I|81(l via the 
quantity 9 — (MV 2 /W Dofi±) . For ablock sliding under 
its own weight (W = Mg), with Dq > 1/xm, fi± « 10~ 2 , 



< lO" 5 (V, 



fim I sec j 



can safely be neglected in our low 
velocity regime, and the position of the bifurcation line 
is simply given by : 



M0 — Mi 



(82) 



For the RR model, M0 — Mi = B~A is a constant, (K/W) c 
is V-independent. 

At the bifurcation, the two (complex conjugate) neu- 
tral modes oscillate at the critical frequency : 



l>q V Mi 



(83) 



That is, the bifurcation is of the Hopf type. A third or- 
der perturbation expansion can be performed standardly 
provided that and /i± are not mere constants (see §11) . 
It shows that the stick-slip bifurcation is direct - i.e. that 
the SS amaplitude grows continuously when decreasing 
e.g. K below its critical value. 



Note finally that, for a velocity-strengthening /i^, (m<^ — 
fi±) < 0, the system is always stable against infinitesimal 
perturbations. This does not preclude t he p ossibility of 
a finite amplitude instability. Brocklcy 108J has shown 
that - in the case where fi^ = - when one takes iner- 
tia into account, the system does exhibit a strongly hys- 
teretic Hopf bifurcation. As explained in §11, one expects 
this situation to prevail, for rough-on-rough MCI, at large 
enough driving velocities (V > V m i n ) for which geomet- 
ric aging becomes ina ctive . A Brockley-like regime has 
indeed been observed |l09j |. for a paper/paper interface, 
in the cm/sec range. 



Appendix D : The Schallamach model of 
adsorption-controlled friction 

Following Schallamach's seminal article :64], consider 
an extended interface between a soft slider and a hard flat 
track covered with adhesive sites which can pin the slider 
molecules. Represent the slider as a set of Nq identical 
and independent chains of stiffness n, potentially form- 
ing bonds by adsorption of their end monomer onto the 
track. Adsorption and subsequent desorption are ther- 
mally activated, and, when sliding, desorption is aided 
by advection. The number N of bonds is therefore ex- 
pected to be a function of the sliding velocity V and of 
the temperature T. In steady motion, the elastic force 
exerted on a given bond increases linearly with time until 
the bond snaps off and the stored energy is dissipated. 

The frictional force thus reads : 



F = NnVt 



(84) 



where i is the average lifetime of a bond. 

Let r be the average time for which a chain remains 
depinncd before readsorbing. The stationary number of 
bonds is : 



N = N 



and the friction force 



F = NokVt 



t + T 



(t/l 



(85) 



l + (i/r) 



(86) 
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FIG. 34: Sketch of the adhesive pinning potential (full line). 
When the pinned chain is stretched, the desorption barrier 
decreases from AE = W + E to AE* (dashed line) . 



In order to desorb (resp. adsorb), a chain must over- 
come an energy barrier of height AE — W + E (resp. 
E) , as sketched on Figure When the pinned chain 
is stretched at velocity V, the barrier is lowered and its 
effective height becomes AE* — AE — nbVt, where b is 
a length of atomic order. Schallamach assumes that the 
resulting desorption rate is : 



r(t) = — exp 

TO 



AE*(t) 



knT 



(87) 



where Tq -1 is an attempt frequency, and t the age of the 
bond. 

Note that this expression for r tacitly assumes that 
thermal activation is efficient enough for all bonds to 
break before reaching the deterministic threshold such 
that AE* = 0, i.e. that AE*(t) > k B T. We will discuss 
this assumption in more detail below. 

The time r is that for thermal activation above the 
barrier E, of unspecified origin, and Schallamach sets : 



where r out = To exp (AE /k B T) is the desorption time of 
the unstretched chain, and Vq = k B T '/ ' Kbr OU f 

This, together with expression l|86(l . yields the friction 
force in the steady state. More precisely, in the small and 
large velocity limits : 



t/r 



i/r - V ° 

1 1 ~out — "T7 



V_ 

V> 



V 

In — -0.577 - 
Vo 



(V « Vo) 



(91) 



(V » Vo) (92) 



t decreases monotonously with V, while the elastic en- 
ergy nbVt stored before debonding increases from linearly 
(V <§; Vo) to logarithmically (V ^> Vo) as shown on Fig- 
ure 1351 Also shown on this figure is the average number 
of bonds in steady state N(V), which decreases all the 
more slowly to zero in the large V limit that r/r OM t is 
small. 




t = r exp (E/k B T) 



(88) 



FIG. 35: Qualitative plot of the average pinning force f p (thin 
line), the fraction of bonded chains N/No (dashed line), the 
friction force ^(V) (thick line) as predicted by Schallamach's 
model, for r/T out < 1. 

As a result F(V) exhibits a maximum at V = V max 
(Figure |2SJ resulting from the interplay of two effects : 

- the decrease with V of the number of bonds. 

- the increase of the average pinning force f p = nVt. 
In the limiting regimes : 



an expression which neglects in particular advection ef- 
fects. 

Consider a set of no bonds all formed at the same time 
t = 0. At time t, the number of these, n(t), which have 
not desorbed obeys : 



dn 

M = - r{t)n 



(89) 



The average bond lifetime t then reads, with the help 
of equations |8T )l .l|%9" )l : 



t= / dt = T out — / — , nr /T » dy 90 

Jo n o V J y+ (Vo/V) 



F(V) ~ N oK - 



-V 



Tout 



y ' t b V Vo 



(V « Vb) (93) 



{V->V max ) (94) 



Let us now discuss in more detail Schallamach's results 
and assumptions. Two weak points in the theory are 
immediately clear. 

• It predicts vanishing friction at large V, since, in 
this limit, complete depinning is achieved. The junc- 
tion is then a mere sheared liquid layer, whose viscous 
dissipation has been neglected. The corresponding con- 
tribution F v is(V) to the total friction force F to t(V) = 



3G 



F(V) + F vis (V) becomes dominant for V ^> V max . Its 
precise expression of course depends on the nature of the 
slider material (e.g. hydrogel vs elastomer). 

• As mentioned above, Schallamach's expression for 
the desorption rate r is valid only as long as AE*(t) <C 
ksT. Since equation (|74J) results in an unboundedly 
growing value of nbVt (see Figure I35fl , this assumption 
fails for V > V* such that AE = nbV*t(V*), which 
yields V* ~ Vq exp(Ai?/fcsT). In this fast regime, we 
are back to the scenario of § II. C. 4. Advection-controlled 
deterministic bond breaking becomes more efficient than 
"premature" thermally activated depinning, resulting in 
a value of f p which crosses over from the Schallamach 
logarithmic regime to the deterministic saturation value 
fsat = AE/nb. That is, as shown on Figure 03 f P ex- 
hibits three regimes : an initial linear one for V < V, 
a logarithmic one for V < V < V* , saturation for 

V > V*. The lower crossover V can be evaluated as 

V =s V {AE/k B T) ~ AE/nbrout. 

Keeping track of these two corrections results in 
a total friction force curve whose qualitative shape 
depends on the order of magnitude of the parameter 
t I Taut. Note that, within Schallamach's assumption 
(equation T/r out < 1. 



If t I T out < 1, one easily evaluates V max w V Q T out /T, 
so that 



V 



AE 



Vrr, 



k B T r„ 



AE 



-W/k B T 



< 1 



and F(V) exhibits a wide plateau followed by a slow de- 
crease. Whether or not the total force F to t = F + F V i S 
exhibits a A/"-shape for V > V max depends on the effec- 
tive viscosity of the unpinned junction and on the density 
of pinning sites N . 

As t I T ou t increases towards 1, the width of the plateau 
decreases with V max , F(V) becomes steeper on the high- 
ly side, so that the more likely a Af-shaped F tot {V), hence 
a ^-weakening friction regime leading to stick-slip. 

Finally, it is worth recalling that, when the system 
is loaded from rest at a prescribed velocity Vi oa d larger 
than, typically, V, its response is initially quasi-elastic, 
the upper limit of this regime appearing as a static 
threshold. Since V ~ AE / K.bT ou t, reasonable estimates 
lead to conclude that in order for V to lie in the sub- 
/j,m/sec range, adsorption should be quite strong - cor- 
responding to binding energies W not far below the eV 
level. 
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